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1 .  INTRODUCTION 


Beginning  with  the  simultaneous  discovery  of  the  fundamental  algorithm 
by  Hirota,  et  al.  (1970),  and  by  Roache  (1971),  a  new  class  of  algorithms  for 
solving  elliptic  partial  difference  equations  has  been  developed  during  the 
last  10  years.  Roache  (1971)  noted  that  the  fundamental  algorithm,  which  is 
analogous  to  "shooting  methods"  in  ordinary  differential  equations,  and 
described  as  the  "error  vector  propagation"  (EVP)  method,  is  very  efficient 
when  applied  to  low-resolution  problems,  but  cannot  be  applied  to  high- 
resolution  problems  due  to  roundoff  error.  Dietrich  (1974)  discovered  an 
iterative  application  of  the  algorithm  that  could  be  applied  to  high- 
resolution  problems,  and  this  iterative  technique  was  further  explored  and 
refined  by  Dietrich,  et  al.  (1975)  and  by  Dietrich  (1977).  Madala  (1978) 
discovered  a  method,  called  the  "stabilized  error  vector  propagation"  (SEVP) 
technique,  whereby  the  basic  algorithm  can  be  applied  as  a  direct  method  for 
high-resolution  problems.  Roache  (1978)  reviewed  the  development  of  EVP  and 
analyzed  its  characteristics  in  detail. 

One  important  class  of  applications  of  EVP-based  elliptic  solvers  is 
in  semi-implicit  formulations  of  the  fluid  dynamic  equations  in  which  dominant 


I.  Hirota,  T.  Tokioka,  and  M.  Nishiguchi,  J.  Meteor.  Soc.  Japan  ji8^  (1970), 
pp.  161-167 . 

P.  J.  Roache,  "A  New  Direct  Method  for  the  Discretized  Poisson  Equation,"  in 
Proceedings  of  the  Second  International  Conference  of  Numerical  Methods  in 
Fluid  Dynamics,  Berkeley,  California,  September  16-19,  1970,  pp.  48-53,  1971. 

D.  Dietrich,  J.  Meteor.  Soc.  Japan  ^2^  (1974),  pp.  337-338. 
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pp.  421-439. 

D.  Dietrich,  Final  Report,  Contract  No.  N00014-77-C-0208  for  the  Naval 
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R.  Madala,  Mon.  Wea.  Rev.  106  (1978),  pp.  1735-1741. 

P.  J.  Roache,  Nura.  Heat  Trans.  1  (1978),  pp.  1-201. 


terms  involved  in  gravity  (O'Brien  and  Hurlburt,  1972)  and/or  compression  wave 
(Dietrich,  et  al.,  1981)  propagation  are  solved  implicitly  (with  the  remaining 
terms  explicit) .  A  more  powerful  application  exists  in  solving  "almost  fully 
implicit”  formulations  of  the  fluid  dynamic  equations  (Dietrich,  1975);  this 
application  is  recommended  when  the  shortest  fluid  dynamic  time  scale  of 
interest  is  long,  compared  to  the  well-known  CFL  time  step  limit,  as  appears 
to  be  the  case  in  most  fluid  dynamic  problems  with  constant  or  slowly  varying 
external  forcing.  (In  some  problems  with  waves  generated  by  time-varying 
external  conditions,  this  might  not  be  the  case.) 

In  Section  2  of  this  report,  we  describe  an  SEVP  algorithm,  program 
JAYNOR,  for  solving  general  two-dimensional,  second-order  elliptic  equations. 
This  algorithm  can  be  applied  to  linear  problems  with  variable  coefficients, 
irregular  boundaries,  islands,  and  general  boundary  conditions.  Its  limita¬ 
tions  are  mainly  associated  with  roundoff  error,  as  described  in  Section  3. 
In  Section  4,  we  describe  how  to  use  a  computer  program  based  on  this  general 
algorithm,  as  well  as  more  specialized  versions  of  this  program.  Appendix  II 
contains  listings  of  this  general  program  and  its  more  specialized  versions. 


J.  J.  O'Brien  and  H.  E.  Hurlburt,  J.  Phys .  Oceanogr .  (1972) ,  pp.  14-26. 

D.  Dietrich,  H.  Klein,  and  R.  K.-C.  Chan,  "FLAME:  A  Computer  Model  of  Time- 
Dependent,  Multidimensional,  Multiphase  Reactive  Flow,"  paper  presented  at  the 
20th  National  Heat  Transfer  Conference,  Milwaukee,  Wisconsin,  August  2-5, 
1981. 


D.  Dietrich,  J.  Meteor.  Soc .  Japan  53  (1975),  pp.  222-225 


2.  THE  GENERAL  SEVP  ALGORITHM 


The  general  SEVP  algorithm  described  here  Is  similar  to  the  one 
described  by  Madala  (1978).  The  differences  are  mainly  the  capability  of 
Including  Islands  In  the  present  version,  and  the  elimination  of  the 
homogeneous  solution  calculation  during  the  forward  sweep  by  more  fully 
utilizing  influence  coefficient  arrays  generated  during  a  preprocessor. 

The  new  SEVP  algorithm,  JAYNOR,  uses  a  master  flag  array,  N^ 
defined  by: 


N 


1  if  the  (i,j)  point  is  "inside  the  boundary" 

0  otherwise 


Points  on  islands  are  considered  "outside  the  boundary".  If  Nj  j  =  1,  the 
following  equation  is  satisfied: 


A  a  *  X.  . +  B,  *  X.  .  _  ...  +  C.  •  X...  . 

i,j  i,j+l  i,j  i+2,j+l  i,j  i+l,j 


+  Di,j  *  Xi+l,j+2  +  Ei,j  *  Xi+l,j+l  “  Fi,j  » 


(1) 


where  Aj  j,  j,  j,  Dj  j,  and  E^  j  are  given  coefficient  arrays  of  the 
five-point  elliptic  operator;  Fj  j  is  a  given  source  term  array;  and  j  is 
the  solution  array  determined  by  JAYNOR.  If  N^  j  ■  0,  the  j  is  left 
unchanged.  In  addition  to  the  above-indicated  parameters  and  others  normally 
needed  to  define  a  two-dimensional  elliptic  partial  difference  equation  with 
second-order  accuracy,  JAYNOR  requires  the  user  to  specify  the  number  of 
subregions  in  which  the  aforementioned  SEVP  algorithm  is  to  applied,  and  the 
grid  lines  separating  these  subregions.  In  general,  the  more  subregions  used, 
the  smaller  the  roundoff  error  will  be  in  the  solution  (at  the  expense  of  more 
auxiliary  storage  and  calculation).  These  grid  lines  are  transverse  to  the 
EVP  march  direction,  which  is  the  j-direction  in  JAYNOR. 

We  now  describe  the  JAYNOR  algorithm  in  detail.  Since  Neuman  and 
mixed  boundary  conditions  can  be  easily  transformed  into  Dirichlet-type 
boundary  conditions  by  suitably  redefining  the  partial  difference  operator  and 
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source  term  at  points  adjacent  to  boundaries  (see  Appendix  I),  we  describe  its 
application  to  problems  With  Dirichlet  conditions.  (JAYNOR  assumes  the  user 
has  already  performed  this  simple  task.) 

First,  we  describe  JAYNOR  application  to  a  regular  region, illustrated 
in  Figure  1.  Two  classes  of  auxiliary  influence  arrays  are  generated  in  the 
JAYNOR  preprocessor,  subroutine  JAY.  Both  classes  are  composed  of  selected 
subvectors  from  the  total  solution  vector,  forced  by  a  selected  isolated  unit 
source  term.  In  a  regular  problem,  each  such  subvector  of  the  first  class  is 
the  solution  at  the  second  row  from  the  bottom  of  a  subregion,  forced  by  a 
unit  source  term  assigned  to  a  selected  point  in  the  second  row  from  the  top 
of  the  same  subregion,  in  an  associated  problem  in  the  region  below  the  top 
boundary  of  the  subregion  —  the  associated  problem  has  homogeneous  boundary 
conditions,  including  at  the  top  bounday  of  the  subregion,  and  has  zero  source 
term  everywhere,  except  at  the  selected  unit  source  point.  In  JAYNOR,  this 
array  is  called  RINV;  specifically,  RINV(I1 ,I2,K)  is  the  solution  at  the 
(12+1)  position  of  the  second  row  from  the  bottom  of  subregion  K,  forced  by  an 
isolated  unit  source  at  the  (11+1)  position  in  the  second  row  from  the  top  of 
subregion  K.  In  a  regular  problem,  each  subvector  of  the  second  class  Is  the 
solution  at  the  second  row  from  the  top  of  a  subregion,  forced  by  a  unit 
source  term  assigned  to  a  selected  point  In  the  second  row  from  the  top  of  the 
same  subregion,  again  in  the  associated  problem  described  above.  In  JAYNOR, 
this  array  is  called  RINV1;  specifically,  RINV1(I1,I2,K)  Is  the  solution  at 
the  (12+1)  position  of  the  second  row  from  the  top  of  subregion  K,  forced  by 
an  isolated  unit  source  at  the  (11+1)  position  in  the  second  row  from  the  top 
of  subregion  K.  The  top  subregion  does  not  require  an  array  of  this  second 
class.  (Although  the  above  description  indicates  the  storage  sequence  in 
regular  problems,  arrays  RINV  and  RINVl  are  stored  as  one-dimensional  arrays 
to  save  storage  in  Irregular  problems  due  to  differing  storage  requirements  in 
each  subregion  —  see  below.) 

Using  the  arrays  RINV  and  RINVl,  as  determined  by  the  JAYNOR  pre¬ 
processor,  which  is  performed  in  subroutine  JAY,  subroutine  NOR  calculates  the 
solution  to  the  partial  difference  equation  (1)  with  assigned  boundary 
conditions.  This  is  done  using  a  forward  sweep,  followed  by  a  backward  sweep, 
similar  to  the  procedure  described  by  Madala  (1978),  utilizing  RINV  and  RINVl 
data  along  the  way.  In  the  present  method,  the  forward  "particular  solution 
sweep"  and  backward  "homogeneous  solution  sweep"  each  requires  slightly  fewer 
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Figure  1.  Application  of  program  JAYNOR  to  a  regular  mesh.  The  partial 
difference  equation  is  satisified  at  the  center  of  each  cell, 
except  those  outside  the  boundary,  which  are  shaded;  the  source 
term  and  solution  arrays  are  cell-centered.  Grid  lines  forming  the 
cell  boundaries  are  drawn.  Grid  lines  between  SEVP  subregions  are 
darkened.  In  this  example,  JAYNOR  program  topological  variables 
are  as  follows:  M«14,  N«20,  NSB*3,  IE(1)«7,  IE(2)-13,  IE(3)*19, 
and  the  NTYPE(I,J)  array  has  zeroes  for  1  =  1,  1*14,  J-l  and  J*20, 
and  ones  everywhere  else  (there  are  no  islands).  SEVP  starting 
vector  positions  are  indicated  by  S's,  and  SEVP  residual  vector 
positions  are  indicated  by  R's;  roost  of  the  R's  have  zero  sub¬ 
scripts  to  indicate  they  are  "open  residual  positions"  (see  main 
text) . 
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operations  than  a  single  ordinary  relaxation  sweep,  except  for  operations 
involving  RINV  and  RINV1.  Subroutine  NOR  avoids  generating  a  homogeneous 
solution  on  the  forward  sweep  by  using  RINV1  to  adjust  the  grid  point  values 
in  the  second  row  from  the  top  of  each  subregion  (except  the  top  subregion) 
directly  during  the  forward  sweep,  instead  of  using  RINV,  RINV1,  and  an 
appropriate  homogeneous  sweep  to  create  these  same  adjustments,  as  described 
by  Madala  (1978).  Every  other  step  of  the  present  procedure  is  essentially 
the  same  as  described  by  Madala  when  JAYNOR  is  applied  to  a  regular  region. 

Application  of  JAYNOR  to  elliptic  problems  in  regions  with  irregular 
boundaries  and  islands  is  similar  to  its  application  to  regular  problems.  The 
differences  are  mainly  in  the  definition  and  use  of  the  preprocessor  arrays 
RINV  and  RINVl .  These  differences  are  associated  with  the  definition 
"starting  vector  points”  and  "error  vector  points”  in  each  subregion.  For 
regular  problems,  "starting  vector  points”  are  points  in  the  second  row  from 
the  bottom  of  each  subregion,  and  "error  vector  points"  are  points  in  the 
second  row  from  the  top  of  each  subregion  (see  Figure  1).  When  the  problem  is 
irregular,  "starting  vector  points"  are  all  points  in  the  second  row  from  the 
bottom  of  each  subregion  that  are  inside  the  boundary  (NTYPE=1)  plus  all 
points  inside  the  boundary  that  are  immediately  above  points  that  are  outside 
the  boundary  (NTYPE=0)  (see  Figure  2).  Similarly,  "error  vector  points”  are 
all  points  In  the  second  row  from  the  top  of  each  subregion  that  are  inside 
the  boundary  plus  all  points  inside  the  boundary  that  are  immediately  below 
points  that  are  outside  the  boundary  (see  Figure  2).  The  general  defintions 
for  RINV  and  RINVl,  applicable  for  both  regular  and  Irregular  problems,  are 
then  as  follows.  RINV(L,M,K)  is  the  solution  at  the  Mth  starting  vector 
position  in  subregion  K,  forced  by  an  isolated  unit  source  at  the  error 

vector  position  in  subregion  K,  in  an  associated  problem,  as  described 
above.  RINV1(L,M,K)  is  the  solution  at  the  "open"  error  vector  position 

in  subregion  K,  forced  by  an  isolated  unit  source  at  the  error  vector 

position  in  subregion  K,  in  the  associated  problem,  where  "open"  means  that 
the  point  above  the  error  vector  position  is  inside  the  boundary  (NTYPE=1) 
(see  Figure  2).  Since  the  ranges  of  L  and  M  can  vary  from  subregion  to 
subregion  in  irregular  problems,  arrays  RINV  and  RINVl  are  treated  as  one¬ 
dimensional  arrays  to  save  storage. 
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Figure  2.  Application  of  program  JAYNOR  to  an  irregular  mesh  with  islands. 

Refer  to  Figure  1  caption  to  interpret  this  figure.  Some  of  the 
R's  have  zero  subscripts  to  indicate  these  are  "open  residual 
positions"  (see  main  text).  The  user  need  not  specify  the  R,  S, 
and  Rq  positions;  they  are  determined  in  subroutine  TOPOL  from  the 
user-specified  topological  flag  array  NTYPE  (see  main  text). 


3.  SELECTION  OF  MARCHING  DIRECTION  AND  BLOCK  SIZE 
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The  choice  of  marching  direction  and  block  size  is  determined  by  the 
accuracy  desired  in  the  solution,  the  computing  precision  used,  and  the  EVP 
roundoff  error  characteristics  of  the  partial  difference  operator.  (In  many 
applications,  it  is  desirable  to  use  double  precision  on  IBM  computers  and  on 
other  computers  for  which  double-precision  computation  is  about  as  fast  as 
single-precision  computation.  This  can  usually  be  facilitated  by  using 
IMPLICIT  REAL*8  type  statements  at  the  beginning  of  the  EVP  subroutines.) 

Detailed  discussions  on  this  are  given  by  Dietrich  (1975)  and  Roache  (1978). 
In  general,  the  choice  of  march  direction  and  block  size  depends  on  a  number 
of  factors,  and  only  rough  guidelines  can  be  given.  Fortunately,  this  is  all 
that  is  needed,  since  even  when  operating  well  within  roundoff  error  restric¬ 
tions,  the  SEVP  method  is  robust  and  is  usually  competitive  with  the  best  of 
other  methods,  often  being  superior.  However,  it  should  be  noted  that  when 
the  geometry  is  regular  (or  nearly  regular),  and  the  partial  difference 
operator  has  constant  coefficients,  there  are  transform  methods  (combined  with 
a  capacitance  matrix  to  handle  any  irregularities)  that  can  be  superior  to 
SEVP  in  terms  of  storage  requirements  and  competitive  in  terms  of  operation 
count  (Wallcraft,  1980). 

A  general  rule-of-thumb  guideline  for  the  choice  of  subregion  size 

when  using  SEVP  is  that  for  an  elliptic  partial  differential  equation  with  low 
to  moderate  diagonal  dominance  (as  in  a  Poisson  partial  differential 
equation),  the  normalized  absolute  error  In  the  application  of  the  EVP 

T> 

algorithm  to  a  given  subregion  is  of  order  5  times  the  computing  precision 
used,  where  B  Is  the  ratio  of  the  subregion  length  in  the  EVP  march  direction 
to  the  minimum  grid  interval  size  in  the  other  direction.  For  example,  for  a 
Poisson  equation  on  a  uniform  mesh  (Ax  *  Ay),  and  computing  precision  of  15 
significant  figures,  the  SEVP  solution  would  be  accurate  to  about  five 

significant  figures  if  the  maximum  subregion  length  were  about  14  grid 
Intervals.  However,  SEVP  can  be  applied  iteratively  (as  can  any  direct 

solver)  to  reduce  roundoff  error  to  any  desired  level  above  the  computing 


A.  Wallcraft,  "Capacity  Matrix  Techniques,"  Ph.D.  Thesis  (1980),  Department  of 
Mathematics,  Imperial  College,  London,  England  S.W.7. 
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precision  available.  For  example,  if  the  roundoff  error  is  0(10”'*)  after  one 
call  to  SEVP,  the  error  can  be  reduced  to  0(10“^®)  after  two  SEVP  calls.  Each 
succesive  SEVP  problem  is  forced  by  the  roundoff  error  residuals  from  the 
previous  SEVP  problem,  and  the  result  is  added  to  the  previous  solution. 

A  rule  for  choosing  the  SEVP  march  direction  is  suggested  by  Madala 
(1978)  .  Madala  recommends  choosing  the  direction  for  which  "A  >  1",  where  A 
is  the  ratio  of  the  grid  interval  transverse  to  the  SEVP  march  direction 
divided  by  the  grid  interval  in  the  SEVP  march  direction.  However,  in  some 
problems  there  are  regions  for  which  A  >  1  and  other  regions  for  which  A  <  1. 
In  such  problems  the  use  of  SEVP  as  a  direct  method  may  be  costly.  A  good 
possible  alternative  in  such  cases  is  to  use  SEVP  iteratively  (with  variable 
march  direction)  in  the  same  way  EVP  was  used  iteratively  by  Dietrich,  et  al. 
(1975);  in  some  cases,  it  may  be  useful  to  use  alternate  forms  of  EVP  such  as 
described  by  Dietrich  (1977)  and  by  Roache  (1978).  Even  in  problems  for  which 
"A  -  1“  has  the  same  sign  everywhere,  it  might  not  be  best  to  march  in  the 
direction  for  which  A  >  1.  Specifically,  when  A  >  1  for  a  march  direction 
with  much  fewer  mesh  points  than  in  the  other  direction,  the  other  march 
direction  may  be  better  even  though  more  SEVP  subregions  are  required,  since 
the  influence  arrays  RINV  and  RINV1  will  be  much  smaller  for  a  given 
subregion.  This,  of  course,  can  be  affected  by  the  kind  of  computer  being 
used  (e.g.,  vector  vs.  scalar  computer). 

Finally,  when  there  are  islands,  it  is  sometimes  advantageous  to  place 
the  SEVP  subregion  boundaries  in  such  a  manner  that  the  islands  straddle  the 
boundaries,  especially  when  this  does  not  increase  the  number  of  SEVP 
subregions  required;  this  reduces  the  number  of  starting  and  error  vector 
points  required  and  the  associated  auxiliary  computation  and  storage  used  by 
the  SEVP  algorithm.  Similarly,  if  there  are  no  islands,  it  is  desirable  to 
place  the  SEVP  subregion  boundaries  near  the  narrowest  parts  of  the  SEVP 
domain,  such  as  the  throat  of  an  hourglass-shaped  domain. 


n 


4.  USING  THE  JAYNOR  COMPUTER  PROGRAM 


Program  JAYNOR  sets  up  and  solves  a  general  sample  problem,  using 
preprocesor  subroutines  TOPOL  and  JAY,  and  repeat  solution  subroutine  NOR;  In 
JAYNOR  program  variables,  the  sample  problem  may  be  written: 

AX(I, J)*X(I, J+l )+CX(I, J)*X(I+2 , J+l )+AY(I, J)*X(I+1 , J) 

+CY(I,J)*X(I+l,J+2)+BB(I,J)*X(I+l,J+l)  =  F(I,J)  .  (2) 

If  NTYPE(I, J)  =  1,  Eq.  (2)  Is  satisfied;  If  NTYPE(I,J)  -  0,  X(I,J)  Is  left 
unchanged  and  treated  as  a  boundary  condition,  where  appropriate.  Thus,  the 
user  must  specify:  the  partial  difference  operator  coefficient  arrays  AX,  CX, 
AY,  CY,  and  BB;  the  source  term  array  F;  the  topological  array  NTYPE;  the  x- 
dimension  of  the  solution  array  M;  the  y-dimenslon  of  the  solution  array  N; 

the  maximum  (for  all  SEVP  block  subregions)  number  of  starting  vector  guess 

positions  MXGS;  the  number  of  SEVP  block  subregions  NBLK;  the  dimensions  of 
the  Influence  coefficient  arrays  (RINV  and  RINV1)  LD1  and  LD2;  and  the  grid 

line  number  at  the  top  of  each  subregion,  array  IE.  If  any  of  the  user- 

specified  dimensions  MXGS,  LD1,  or  LD2  is  too  small  for  the  specified  mesh 
topology  (NTYPE  array),  appropriate  diagnostic  messages  are  printed  in 
subroutine  TOPOL  and  the  computation  Is  halted.  Most  arrays  have  dimensions 
selected  from  the  following  quantities:  M,  M-2,  N,  N-2,  MXGS,  NBLK,  and 

NBLK-1.  Specifically,  the  appropriate  dimensions  are  as  follows: 

(M-2, N-2):  AX,AY,CX,CY,BB,Q,F 
(M,N) :  RECUR, X,H, ERR, XX, NTYPE 
Also:  DUM0( MXGS, NBLK-1) ,DUM1(MXGS) ,DUM2(MXGS) , 

ILFT(N) ,IRGT(N) ,IGES( MXGS, NBLK) ,JGES( MXGS, NBLK) , 

ICHK( MXGS, NB LK) ,JCHK( MXGS, NBLK) ,NEV(NBLK) , 

LINK(M-2 , NBLK-1 ) , NLINK( NBLK-1 ) ,  IE( NBLK) , LI ( NBLK) , 
and  L2( NBLK-1). 
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The  maximum  lengths  required  by  RINV  and  RINV1  are: 

RINV(NBLK*MXGS*MXGS) ,RINVl((NBLK-l)*MXGS*(M-2)) . 

A  complete  alphabetical  list  of  the  JAYNOR  program  significant  variables  and 


their  definitions 

asterisk. 

follows.  User  -specified  variables  are 

TABLE  1:  JAYNOR  VARIABLES  AND  DEFINITIONS 

denoted  by  an 

Variable  Name 

Variable  Definition 

AX(I,J)* 

Left  coefficient  of  difference  operator 
(I+l.J+l) 

at  mesh  point 

AY(I, J)* 

Bottom  coefficient  of  difference  operator 
(I+l.J+l) 

at  mesh  point 

BB(I.J)* 

Central  coefficient  of  difference  operator 
(I+l.J+l) 

at  mesh  point 

CX(I,J)* 

Right  coefficient  of  difference  operator 
(I+l.J+l) 

at  mesh  point 

CY(I,J)* 

Top  coefficient  of  difference  operator 
(I+l.J+l) 

at  mesh  point 

ERR(I,J) 

Normalized  residual  array 

ERSUM 

Mean  absolute  residual 

F(I,J)* 

Source  term  for  difference  operator  at  mesh 

point  (I+l.J+l) 

FSUM 

Mean  absolute  source  term 

H(I,J) 

Alraost-homogeneous  solution  component  array  calculated  on 
backward  SEVP  sweep 

(IC,JC) 

Mid-point  of  island  location  in  SEVP  sample 

problem 

ICHK(N.NB) 

The  N1"*1  residual  position  in  the  NBth  block  is  on  the 
ICHK(N.NB)  vertical  mesh  line 

(continued) 
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TABLE  1  (continued) 


■ - -  - 

Variable  Name  Variable  Definition 


IE(NB)  The  top  azimuthal  mesh  line  inside  the  NBC^  block 

(IE(NBLK)=N-1) 

IGES(N,NB)  The  starting  vector  position  in  the  NBth  block  is  on 

the  IGES(N,NB)  vertical  mesh  line 


ILFT(J) 


The  left-most  point  inside  the  boundary  on  horizontal  mesh 
line  J 


|  IRGT(J) 


JCHK(N,NB) 

JGES(N,NB) 


The  right-most  point  inside  the  boundary  on  horizontal  mesh 
line  J 

The  residual  position  in  the  NBC^  block  is  on  the 

JCHK(N,NB)  horizontal  mesh  line 

The  NT^  starting  vector  position  in  the  NB^  block  is  on 
the  JGES(N,NB)  horizontal  mesh  line 

The  dimension  of  array  RINV 


LD2* 

[ 

LINK(N.NB) 


Ll(NB) 


The  dimension  of  array  RINV1 

The  Nth  open  residual  position  in  the  NBth  block  is  at  the 
LINK(N.NB)  residual  position  ( I=ICHK( LINK(N,NB) ,NB) , 
J-JCHK(LINK(N,NB),NB)«IE(NB)) 

The  total  number  of  elements  in  array  RINV  for  blocks  below 
block  NB 


L2(NB) 


MXGS* 


NBLK* 


The  total  number  of  elements  in  array  RINV1  for  blocks 
below  block  NB 

The  number  of  vertical  mesh  lines  spanning  the  problem  to 
be  solved 

The  maximum  number  of  guess  points  to  be  encountered  over 
all  blocks  (can  be  larger) 

The  number  of  horizontal  mesh  lines  spanning  the  problem  to 
be  solved 

The  number  of  blocks  (SEVP  subregions)  to  be  used 


NEV(NB) 


The  number  of  starting  vector  (or  residual)  positions  in 
block  NB 


(continued) 


TABLE  1  (continued) 


Variable  Name 

Variable  Definition 

NLINK(NB) 

The  number  of  open  residual  positions  in  block  NB 

NTYPE( I, J)* 

Topological  flag  array  ( NTYPE( I, J)*l  if  the  partial 
difference  operator  is  to  be  satisfied  at  the  grid  point 
located  at  the  intersection  of  the  It^1  vertical  mesh  line 
and  horizontal  mesh  line;  NTYPE(I,J)*0  otherwise) 

Q(I,J) 

Same  as  F(I,J) 

RECUR(I.J) 

Real  array  of  l's  and  O' s  related  to  NTYPE 

RINV 

Influence  coefficient  array  relating  unit  source  terms  at 
error  vector  positions  to  solutions  at  starting  vector 
positions 

RINV1 

Influence  coefficient  array  relating  unit  source  terms  at 
error  vector  positions  to  solutions  at  error  vector 
positions 

X(1,J) 

Solution  at  mesh  point  (1,J) 

XX(I,J) 

Same  as  X(1,J) 

I  - 


*User-specif led  variable 
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APPENDIX  I 


TRANSFORMATION  TO  DIRICHLET  BOUNDARY  CONDITIONS  FROM 
NEUMAN  OR  MIXED  BOUNDARY  CONDITIONS 


Using  the  notation  of  Eq.  (1)  in  the  main  text,  suppose  the  point 
(i-I+1, j-J+1)  is  inside  the  boundary,  whereby  Eq.  (1)  applies  with  i=I  and 
j**J,  and  suppose  that  the  point  (i=I,j=J+l)  lies  outside  the  boundary,  whereby 
Eq.  (1)  does  not  apply  with  i"I-l  and  J"J.  Instead,  we  apply  a  "boundary 
condition"  of  the  form 


XI,J+1  “  3  ~  “  XI+1,J+1 


(3) 


where  a  and  3  are  specified.  If  a  ■  0  we  have  a  Dirichlet  condition,  and  no 
adjustment  is  necessary  before  using  program  JAYNOR.  However,  if  a  ^  0  we  can 
use  Eq.  (3)  to  eliminate  Xj  from  the  equation  set.  We  do  this  by 
substituting  Eq.  (3)  into  Eq.  (1)  with  i*I  and  j«J  to  obtain 


AI,J  *  XI,J+1  +  BI,J  *  XI+2 ,  J+l  +  CI,J  ’  XI+1,J 


+  °I,J  ’  XI+l,J+2  +  ^.J  XI+1,J+1  "  *1,3 


(4) 


where  A' 

I,  J 


0, 


a,  and  Fl 

1>  J 


3. 


The  modified  partial  difference  equation  (4)  can  be  solved  by  program 
JAYNOR;  the  value  of  Xj  J+^  does  not  matter  in  this  solution  procedure, 
although  the  user  is  free  to  calculate  the  appropriate  value,  after  JAYNOR 
gives  the  corrct  Xj+^  value,  using  Eq.  (3).  Mixed  boundary  conditions  in 
the  other  three  directions  from  the  point  (I+1,J+1)  can  be  treated  in  an 
analogous  manner. 


APPENDIX  II 


COMPUTER  LISTINGS  OF  FOUR  SEVP  VERSIONS 


The  main  text  describes  the  most  general  version  of  program  JAYNOR  in 
detail.  This  is  Version  2.  The  differences  between  each  successive  version 
and  the  previous  one  are  indicated  by  comment  cards  in  the  main  program 
JAYNOR. 

The  listings  for  four  versions  follow. 


VERSION  2 


1. 

01 

PROGRAM  JAYNORI INPUT .OUTPUT , TAPES=TNPuI ,TAPE6=0UTPUT  1 

2. 

00 

C0MN0N/SEVP/AX<7,9I,AYI7,9I,BB(7,9I,CX(7,9I,CY(7,9I,0(7,9I,F(7, 

1. 

00 

A  RECUR(9,I1I,R1NV(S8BI,RINV1( 1081, 0UMC<I4,2I,0UH1 (141, 

4. 

00 

B  0UM2d4l,X(9,lll,H(9,lll,ERR(9,lII,XX(9,Ill,ILFT(lII, 

S. 

00 

C  IR6T < 11 1 , 16ES <14 ,11 ,JGES <14,11 ,TCHK  < 14, 11 ,JCHK (14 , 11 ,NEV 111 , 

6. 

00 

0  LTNK<7,2I,NLINK<2I,NTYPE<4,1II,TE<5I,LI<1I,L2<2I 

7. 

00 

DATA  N,N,MXGS,NBLK,L0I.L07/9,II,14,1,SBB,108/ 

1. 

00 

OATa  TE/4,7,10/ 

9. 

00 

c 

DIETRICH'S  MODIFIED  HAOALA  SOLVER,  VERSION  2. 

10. 

00 

c 

THIS  VERSION  TS  IDENTICAL  TO  VERSION  1,  EXCEPT 

n. 

00 

c 

FOR  MAKING  RINV  ANO  RINVI  INTO  1 "DIMENSIONAL  ARRAYS,  ANO  MODIFYING 

12. 

00 

c 

BOOKKEEPING  TO  REOUCE  STORAGE  REQUIREMENTS. 

11. 

00 

c 

eouation  to  be  solveo  is 

14. 

00 

c 

NTYPE<I«I.J«II  *  <AX(I,J|  •  X(I,J*II  •  CX(I,JI  •  X(1*2,J*1I 

IS. 

00 

c 

*  AY  (1  ,  J  |  •  Xd'l.JI  *  CY(I,JI  •  X(T«1,J«2II 

14. 

CO 

c 

«  <2  -  NTYPE<I«I,J«1  II  •  BBd.JI  *  X(I*I,J«1I 

17. 

00 

c 

=  2  •  <1  -  NTYPE  <T»I ,  J«1 1 1  *  BBd.JI  *  X<T*1,J*II 

n. 

00 

c 

♦  NTYPE <T«1,J«II  4  FlIiJI,  <<1  =  0,1 . M * 1 1 , ( J=  0, 1 . N- 

l*. 

00 

c 

WHERE  TERMS  OUTSIDE  DIMENSION  BOUNOS  ARE  INTERPRETED  AS  2ER0*S. 

20. 

00 

c 

THUS,  TF  NT  YPE  <  T  ,  J 1  =0 ,  Xd.JI  IS  LEFT  UNCHANGEO. 

21. 

00 

MI=M-I 

22. 

00 

N 1  =N"  I 

21. 

00 

M2=M-2 

24. 

00 

N2=N-2 

2S. 

00 

NBLKlrNBLK-I 

2C. 

00 

00  80  J=1,N 

27. 

00 

00  80  1=1, M 

2o. 

00 

8  C  NTVPE(I,JI=D. 

29. 

00 

00  90  J=2,N1 

10. 

00 

oo  9n  i=2, mi 

11. 

00 

9  C  NTYPE<T,JI=I 

12. 

00 

00  2000  JC=2 ,N I 

11. 

00 

00  2000  1C=2,M1 

14. 

00 

NTYPE  < IC ,JC 1 =0 

IS. 

00 

NTYPE (IC"I , JC 1 =0 

16. 

00 

NTYPE (IC'I , JCI =C 

17. 

00 

NTYPF<IC.JC-1I=0 

18. 

00 

NT  YPE  <  IC,  JC'  1 1=0 

19. 

00 

WRITE <6, 921  < (NTYPE <1, J 1, 1=1, MI,J=I,NI 

40. 

00 

92  F0RHAT(I0X,9III 

41. 

01 

CALL  TOPOL (RECUR, ILFT, IRGT, IGES,J6ES,ICHK,JCHK,NEV, LINK ,NLlNK, 

42. 

00 

A  NTVPE,1E,M,N,MI,N1,H2,NBLK,NBLK1,MXGS,L1,L2,L0I,L02I 

41. 

00 

00  100  J=I,N2 

44. 

00 

00  100  1=1. M2 

45. 

00 

AXd,JI  =  l.*.0I*d"4l**2 

46. 

00 

cx<i,ji=i.-.oi*d-4i4*2 

47. 

00 

AY(I,JI=I.«.0I*<J-5I**2 

48. 

00 

CYd,JI  =  l.-.0I*<J-5|4*2 

49. 

00 

BBd,Ji:"AXd,JI"AYd,JI"CXd,JI-CY(T,JI".I 

SO. 

00 

100  CONTINUE 

SI. 

00 

0EN0»=(.254M4NI**2 

S2. 

00 

00  170  J=I ,N 

SI. 

00 

00  170  1=1, N 

54  . 

00 

Xd,JI  =  d"l,  14  <  I"MI*(  J-I,  |4<  j-ni/oenom 

ss. 

00 

XXd,JI  =  Xd,J|4(l.-NTYPE(T,JII 

56. 

00 

170  CONTINUE 

57  . 

00 

FSUM=0. 
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VERSION  2 
(continued) 


58. 

00 

00  171  J=2,N1 

59. 

00 

L=J-1 

60. 

CO 

00  171  l =2,H I 

61. 

00 

8  =  1-1 

62. 

00 

F|H,LI=AXIH,L19XIX,JI-AYIX,L19XI1,L|.CXI8,LI9XI1.1,J|.CYIX,LI9 

61. 

00 

*  XI 1 , J*1 1 »BB IK.LI9XI 1 , Jl 

6*. 

00 

171 

fsuh=f suh.absifix.lii 

65. 

00 

FSUH=H?6N2/FSUH 

66. 

00 

T  0  =  SE  CONO 1  Tin i 

67. 

01 

CALL  JAYIAX.AY.BB.CX.CY , R 1 N V , R 1 NV 1 , M , RE  CUR , ILF T , 1 RGT , 1 6E 5 , JGE 5 , 

68. 

00 

A1CHX, JCH8.NE V.L1NX .NLINX .NTYPE ,IE ,H . 8 .H2.N2, NBLX .NBLX 1 ,8 XS5 ,L 1 , 

69. 

00 

TlN5UH=0. 

70. 

on 

00  181  J=2,M1 

71. 

00 

00  181  H2.H1 

72. 

00 

011*1. J- 11=F11-1 .  J— 1 1 

71. 

00 

181 

XI1,JI=0. 

7*. 

00 

AL1NE=0. 

75. 

00 

T0=SEC0N0IT1N| 

76. 

01 

CALL  NORIAX.AY.BB.CX.CY.RINY.RINVI.OUMO.OUHI, 0UH2 ,0,M, XX, RECUR, 

77. 

00 

A  lLFT.lRGT.lEES, JGE 5 ,ICKX , JCHX ,NE V ,L 1N8 ,NL1N8 , IE , H,N , H2,N2 .NBLX 

78. 

on 

B  NBLX1,HXGS,L1,L2I 

79. 

00 

DO  990  J=2,N! 

80. 

00 

00  990  1=2. Ml 

tl. 

CO 

9  90 

Xll,JI=Xll.J|.IXIl,JI 

82. 

00 

00  991  J=I,N 

81. 

00 

00  991  1  =  1  .H 

8<t. 

00 

991 

XXI 1 . Jl =0. 

85. 

00 

00  992  J=2,M1 

86. 

00 

L  =  J- 1 

87. 

00 

00  992  1=2. Ml 

88. 

00 

8=1-1 

89. 

00 

9  92 

018, LI =F 18. LI-AYI8.L 19X11, J-1I-AXI8, LI  9X11-1. JI-BBI8, LI  9X11, Jl- 

90. 

00 

A  CX|K,L|9Xll.l,JI-CYIX,LI9Xll,J.]| 

91. 

00 

1C00 

CONTINUE 

92. 

on 

9CS 

F0RHAIIlX.lR15Ct.il 

91. 

00 

rRSuH=o. 

96. 

00 

00  760  1=2, Hi 

95. 

CO 

8  =  1-1 

96. 

00 

00  760  J=2,Nl 

97. 

00 

l  =  J-l 

98. 

on 

ERRI1,JI=FI8,LI-AVIX,LI9XI1,J-1I-AXI8,LI9XI1-1,JI-BBI6,LI9 

99. 

00 

A  XI 1.JI-CXI8.L 19X1 1*1, JI-CY|8,LI*XI1,J«I 1 

ICO. 

00 

ERPll,JI=EBRIl,Jl9NTYREll,JI 

101. 

00 

ERSUH=ER5U»«A65IERRlI.JI 1 

1C2. 

00 

760 

ERRI 1 , Jl =LRRI 1 , JI9F5UH 

101. 

00 

ERSuM=ER5UH9FSUH/|h29N2I 

106. 

00 

WRITE  16,9061 

105. 

00 

9C6 

FORHATI/60X.2CHN0RHALI2E0  PES10UAL5I 

106. 

00 

00  762  lr2, HI 

107. 

00 

WRITE  16,9051 ICBRI 1 , Jl ,J=2,NI 1 

108. 

00 

762 

CONTINUE 

109. 

00 

WR1 TE 1 6 ,906 1  ER5UH 

110. 

00 

9C6 

FORHAT I/J7H  HEAN  AB50LUTE  N0RHAL12E0  RE510UAL  =  .1RE9.2I 

111. 

00 

IF  IER5UH.GT.1.C-6I  5 7 OR 

112. 

00 

IF  ITC.NE.2I  NITRE IlC'l. JCI=1 

111. 

00 

IF  IIC.NE.H1I  NTYPCI1C«1.JCI=1 

116. 

00 

IF  IJC.NE.7l  NTYREI1C,JC-1I=I 

115. 

on 

IF  <JC.NE.N1I  NTYPE  IlC.JCl  1=1 

116. 

00 

2  COO 

NTYRE 1 1C, JCI =1 

117. 

00 

ENO 
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VERSION  2 
(continued) 


1. 

2. 
3. 
6. 

5. 
b  • 
7. 

6. 
9. 

10. 

11. 

12. 
13. 
11. 

15. 

16. 

17. 

18. 

19. 

20. 
21. 
22. 
23. 
29. 

25. 
26* 
27. 

26. 

29. 

30. 

31. 

32. 

33. 
38. 

35. 

36. 

37. 

38. 

39. 
60. 
61. 
62. 
63. 
66. 

65. 

66. 

67. 

68. 
69. 

50. 

51. 

52. 

53. 
56. 

55. 

56. 

57. 

58. 

59. 

60. 
61. 
62. 
63. 
66. 

65. 

66. 

67. 

68. 

69. 

70. 

71. 


02 

00 

00 

00 

00 

00 

on 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

on 

00 

00 

00 

00 

00 

00 

00 

or 

00 

00 

00 

00 

00 

00 

00 

00 

00 

no 

00 

00 

on 

00 

00 

on 

00 

on 

00 

00 

00 

00 

00 

00 


subroutine  jxvixx,xv,B6,cx.cv,rinv,rinvi,h,recur,ilft, 

IIBGT, 16ES, JGCS.1CHK , JCMK ,NE V.LINK ,NLlNK ,NTYPE ,IE , 60  ,N0 ,62 , N2 , 
t  NBLK,NBLKI,"X6S,L1,L2> 

01  PENSION  XXI62.N2) ,XV I62,N2> ,BBf 62.N2I , CX ( 67.N2I , CT 1 62 ,N2) , 

X  R1NVI 1) .RINVI 11 ) , RECUR  I  60 ■ NO) ,6160, NO  I , 

B  ILFTIN0>,1R6T1N0),LIINBLK>,L21NBLKI>,IGESI6XGS,NBLK>, 

C  J6ESC6XGS,N6LK) ,ICMKC6xGS,NBLK I , JCM6  t 6XGS ,NBL* 1 .NEVINBLK ) , 

0  LINK (62.NBLK1 I ,NL1NK (NBLK 1  I ,NT VPE 1 6  0 .NO I , 1 E 1 NBLK I 
JL=I 

LO*L 1 INBLK 1 
NB  =  G 

ICO  NB:N6«1 

NGESXNEVINB) 

NGESHNEVINB-I I 
LXXLl INB I 
LB=L2INB-I I 
IC=L21NB I 
JH=IE INB I 
JMPrjMM 

jhn:jn-2 

00  250  NG:1,N6CS 
1 G  =  IG£  SING ,NB  I 

jg- jge  sing ,nbi 
jf: jg-1 

00  210  JIJL.JHP 
00  210  I -I ,60 
210  N(I,J):0. 

Hlie.JCi:]. 

IF  1JF.GT.JM6)  60  TO  226 
IF  INTVPEilG.JFI.EO. 01  60  TO  220 
N0P=NLINKIN6-1 ) 

00  216  N- 1 ,NOP 
6LL INK (N.NB-l ) 

IF  |ICMK16,NB-1).£0. 16.AN0.JCHK I  6 , NB~I ) .EO.JF )  60  TO  215 
216  CONTINUE 
215  00  21B  N=I,NOP 

LXLINk in.NB-I 1 

2  IB  MIICMXlL,NB-l),jF)=RI*(VllLB*lN-I)*NGESlK6>»CVtI6-I.JG-2) 

220  00  225  JZJF.JM6 

IL-ILFT I J. 2 ) “ 1 
IR  =  IRGTU«?)-I 
222  00  225  1:11, 1R 

225  HII«1,J.2):-RCCURII*1,J*2)*I6X1 1,J)*HI1,J«1).1YII,J)*HII»I,J). 
t  BBII, J)*HII«I ,J«I l«Cxll,J»*H|l<2,J«l ) l/CVIl ,JI 

226  OO  230  NC:l,NGES 
1 : 1CMK INC.NB  )  - 1 
j: JCMK INC.NB )-I 

2  30  R1N»(L»*INC-1)6NGE5«N6i:*XI1 , J I *6 1 1 , J*1 1 *XY 1 1 , J)*m 1 1 • I . J> *BB 1 1 , J» • 

6  HI I«I , J*I l»CXI I , JI»H|I .2, J»1 ) 

IF  INB. EO. NBLK)  60  TO  250 

NOP=NLINKINBI 

jr  IE  I  NR  I 

00  260  NXI.NOP 

6IL INK IN.NB I 

260  R1NYIL0»IN-I I•NGE5♦N6>:MI1CMK(6,NB),J) 

2  30  HUG,  JG  1=0, 

C6LL  6»T1NV|RINVIL»*1I,NGE5,N6ESI 
IF  INB. EO. NBLK)  RETURN 
00  260  i:i,N6E5 
00  260  J=I,NOP 
RINVltLC'l J-1>*NGES«1>=0. 

00  260  K=I,N6E5 

2 60  RINV1 ILC«I J-l 1*N6E5»I ) IR1NY 1 1 1C* I J-l )*NGE5*1 ) - 
6  RINVILXMK-1  )»NGES.I  )*RINYILO«l  J-l  >*NGES*KI 
C  BINVIII.J.NB)  IS  THE  ’  ’  6L60SI  HOMOGENEOUS”  SOLUTION  XT  THE  J-TM  OPEN 
C  RES10U1L  POSITION  FORCEO  BY  X  RESIOUXL  VXLUE  OF  1  XT  The  1-TM 
C  RESIOUXL  POSITION,  MONOGENEOUS  B.C.»5  XRF  X55U6£0  E  VERVwHERE , 

C  INCLUDING  THE  TOP  OF  THE  PRESENT  SUBRE610N  NB . 

JL=JM 
GO  TO  100 
ENO 


22 


VERSION  2 
(continued) 


1. 

01 

SUBROUTINE  NORIAX.AY.BB.CX ,  C Y , B IN  V , R 1 N Y 1 , OUM 0 , 0U« 1 , 0UN2 ,F , N , X , 

2. 

00 

1RCCUR,1LFT,1BGT,16CS,JG£S,1CHK,JCHK,N£V,LINK,NLINK,1C,mo,NO,N2, 

J. 

00 

6  N2,NBLK,N6LK1,MX6S,L.,L?> 

9. 

00 

DIMENSION  1IIN2,N2>,1YIN2,N2>,BBIH2,N2>.CIIH?,N2>,CYIH2,N2>, 

S. 

00 

A  RINVIlI.RINVIll I.RECURlHO.NOl.MIMO.NOI, 

6. 

00 

B  ILFTIN0),1RGTIN0),L1IN8LK)  .L2INBLK 1 I ,1GES INXCS.NBLK 1 . 

T. 

00 

C  JGES INXCS.NBLK )  .1CHK INXGS.NBLK) , JCHK 1 NX6S ,NBLK ) .NEVINBLK > , 

1. 

00 

0  LINKIM2.N6LK1 ) .NLINKINBLKI 1 .IClNBLKt 

9. 

00 

DIMENSION  OUNQ 1 H2, NBLK  H,0UNl|MxGS),0UN2INXGSI,r  IN2.N2l.il  NO.  NO) 

10. 

00 

C0NN0N/BSN/A1T. ALINE. BIT, 8L1NE.1TNAX 

11. 

00 

JS  =  1 

12. 

00 

00  150  NB=1,N6LK 

11. 

00 

LC-L2 1 NB ) 

19. 

00 

JF:I£INB>-2 

IS. 

00 

00  IDS  JXJS.JF 

16. 

00 

IL=1LFT| J«2)-l 

IT. 

oo 

IR-1RBT  I J«2)*l 

16. 

00 

00  IC5  IrU.lR 

19. 

00 

1C5 

XtI.l,J.2)iRECURIIM,J*2>»IFII,J)-AXII,J)»Xll,J«l)-AYIl,J)» 

20. 

00 

A  XI1*1.J>-BBI1,J)*XI1«1,J*1)-CXI1,J1*XII*2.JMII/CYI1,JI«I1.- 

21. 

00 

B  RCCURII«1,J*2))*XI1*1,J*2I 

22. 

00 

IF  IN6.E0.NBLKI  60  TO  150 

21. 

OP 

ngCS=ne  v i nb ) 

29. 

00 

00  Its  n: 1 ,nge s 

25. 

00 

imchkin.nbi-i 

26. 

00 

jrjCHKIN.NItl'l 

21. 

00 

1  IS 

OuNllNI:riI,J)-AXII,J)»Ill,J.l)-lriI,J)»xll.l,J)-BBIl,J)* 

26. 

00 

A  Itl4  1,J4l)-CXIl,J)»ItW,.!«l)-CYtl,J>9XIIM,J«2l 

29. 

00 

NOP-NLlNKINB) 

1C. 

00 

JllEINB) 

11. 

00 

00  120  N=1,N0F 

12. 

00 

0UH2INI=0. 

11. 

00 

00  116  Nrl ,N6ES 

19. 

00 

1  16 

0UH2INl:0UH2  INXOUHl  IN>**1NV1ILC«IN-1>*NGCS«M) 

15. 

00 

HSLINKIN.NB) 

16. 

00 

IzICHKIN.NB) 

IT. 

00 

OUNO 1 N ,NB ) -  X 1 I . J ) 

1*. 

00 

120 

Xll.JI-Xll.JI -0UN2 IN ) 

19. 

00 

150 

JS -IE  1 NB ) 

90. 

00 

00  IPO  NBirl.NBLK 

91. 

00 

nb=nblk-nri»i 

92. 

00 

JSri 

91. 

00 

IF  IN8.NE.1I  JS=1EINB-1> 

99. 

00 

JFSIE 1 NB 1 *2 

95. 

00 

L  A-L 1 1 NB ) 

96. 

00 

L6-L2I NB-1 ) 

9T. 

00 

NGES-NCVINB) 

96. 

00 

IF  INB.EO.NBLKI  GO  TO  201 

99. 

OP 

j:  IEI NB 1 

5G. 

00 

nop  =  nl Ink  I nb ) 

51. 

00 

00  230  N - 1 , NOF 

52. 

00 

nxlinkin.nbi 

51. 

00 

hichkih.nbi 

59. 

00 

2  CO 

xn.jnoUNCiN.NBi 

55. 

00 

2  Cl 

CONTINUE 

56. 

00 

N:1E 1 NB 1 

5T. 

00 

00  202  JTJ5.N 
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VERSION  2 
(continued) 


58. 

00 

00  202  I - 1  •  HO 

59. 

00 

202 

Mf  I v J)=0. 

60. 

00 

00  210  N-l ,NGES 

61. 

00 

I=ICHK IN.NBl-l 

62. 

00 

JXJCHk | N  »NP l-l 

6  3. 

00 

2  10 

OUHlfhIZFf I , JI-AXI I,J>*X(I ,J*1 l-AYf I I«l, JI-BBII, J>* 

64. 

00 

A  X  1 1 ♦ 1  ,JM l-CXI I fJI«X(I«2i JMI-CVII *J>*XII«1*J«2> 

6S* 

00 

00  2 20  N-l *NGES 

66  . 

00 

0UM2I NISO. 

67. 

00 

00  218  M-l ,  NGES 

68. 

00 

2  IS 

0U»*2f  N  >=0UM2  (NMOuMlf  H  >4PIN\f  ILA«IN-1  >4NGES«H> 

69. 

00 

I SIGE  S  f  N ,NP ) 

T  C  • 

00 

J=JGf SfN,NB> 

Tl. 

00 

HII, J|=0UM2I Nl 

T2. 

00 

2  20 

X(I9JUX(I,JI«0UH2fNI 

T3. 

00 

IF  (NB.EQ.l)  GO  TO  2S0 

T4. 

00 

00  22 2  N=1,HX6S 

TS. 

00 

2  22 

0U4 1 f  N 1 -0. 

T6. 

00 

N0P-4L INK INB“1 1 

77. 

00 

J=JS 

T8. 

00 

00  230  N-ltNOP 

79. 

CO 

Hri  INMN.NB-ll 

AC. 

00 

I=ICHKfH,NP-l> 

81. 

00 

2  ?0 

OUHHHI  =  HII,J*l  l*CYl  1-1  ,J-Il 

82. 

00 

C  WHEN  BOUNDARY  IS  IRREGULAR,  THE»E  IS  SOHE  WASTEO  CALCULATION  HERE 

83. 

00 

C  OROER  TO  AxLIO  LOGICAL  QECISIONS  OR  EXTRA  STORAGE  IN6ES  COULO  BE 

64. 

00 

C  REPLACEQ  BY  NLINKINBII 

9S. 

00 

N6T  S  =  NEV INB-1 1 

86. 

00 

00  2A0  N=1,N0P 

BT. 

00 

0U“2INI=P. 

88. 

00 

00  2IS  H=1 , NGES 

89. 

00 

2  IS 

0UH2INI=0uH2IN)«0UHltM>*RINVl(LB*IN>ll*NGES*H» 

90. 

00 

HRL INK  IN ,NB-1 1 

91. 

oo 

2  AO 

HI ICHKI M.NB-l I , JCHK I  M.NB-l  >  >:0UH2IN> 

92. 

00 

210 

00  ICO  jrjS.JF 

93. 

on 

ILIILFTI J»2I-1 

94. 

00 

IRTIRGYI J«2» -1 

9S. 

00 

00  IPO  irlL.lR 

96. 

00 

0UMim:-RECURII*l,J»2>*IAXII,J>*HII,J»l>*AYII,J>*HII»l,JI« 

97. 

00 

A  BBII,JIAHII*1,J»1)»CXII,J)AH(I»2,J»1U/CY(I,J> 

98. 

00 

HIIM,J*2»X0UH1II»*«1.-RECURIIM,J*2M»HIIM,J*2» 

99. 

00 

SCO 

XII*l,J»2irxiI*ItJ»2»«0U«llI » 

00. 

00 

ENO 
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VERSION  2 
(continued) 


1. 

Cl 

SUBROUTINE  TOPOLI RECUR ,ILFT ,IRG1 , ICES , JGES . 1 CHK , JCHK ,NE V .LINK , 

2. 

00 

1Nl1NK,NTTPE,IE.H,N,N1,N1,M2,NBLK,NBLK1,NXGS,L1,l2,L01,L021 

J. 

00 

DIMENSION  RECURIH.NI.  ILF! INI .1RGT INI, IGEStMXGS.NBLK) . 

«. 

00 

A  jGEStHXGS.NBLKI.lCHKlHAGS.NBLK 1 , jCHK < MXGS .NBLK 1 ,NCV (NBLK 1 , 

S. 

00 

BLINK  1 M2.NBLK 1 ),NL INK  1 NBLK 1 1, N T YPE tM,N), IE t NBLK 1,11 1 NBLK ),L2lHBLKl> 

4. 

00 

00  10  J=2,N1 

7. 

00 

HO 

8. 

00 

2 

1=1*1 

9. 

00 

IF  INTYPEtI.J)*NTtPC«l,J-Il.rO.O.ANO.I.LT.Hl)  60  10  2 

10. 

00 

ilftiji=i 

11. 

00 

1  =  « 

12. 

00 

9 

1  =  1-1 

11. 

00 

IF  INTYPEII«JI9NTYPEI1«J-1).EQ.O.ANO.I.GT.I)  60  TO  9 

19. 

00 

1C 

1R6TI J>=I 

IS. 

00 

1R61I 11=2 

16. 

00 

ilfiii>=2 

IT. 

00 

ilftin>=2 

16. 

00 

1r6TIN>:2 

19. 

on 

00  IS  J=2.N1 

20. 

00 

IF  IlLFTIJl.LE.lRSTIJ))  GO  10  IS 

21. 

00 

1RGTI J»=2 

22. 

00 

lLFTI J»=2 

21. 

00 

IS 

continue 

29. 

00 

JM1N=2 

2$. 

00 

HXLNK=0 

26* 

00 

llll>=0 

2T. 

00 

L2< 11=0 

28. 

00 

00  ICO  NB= 1.N8LK 

29. 

00 

NGES=0 

10. 

00 

NCHK=0 

11. 

00 

JHAxriEINB) 

12. 

00 

00  90  1=2, Ml 

11. 

00 

NH=NTtPE 1! .JH1N-1 t 

19. 

00 

NU=NTTPE I  1 ,JH1N 1 

IS. 

00 

00  90  J= JHlN, JNAX 

16. 

00 

NL=NM 

IT. 

00 

NH=NU 

18. 

00 

RECUR  II,J)=NN 

19. 

00 

IF  tNL.E0.01  RECUR 1 1 , J 1=0. 

•a. 

oo 

NU=NHPrU,J*ll 

91. 

00 

IF  INH.EO.Ol  60  TO  90 

92. 

00 

IF  INL.EO.l.ANO.J.NE.JMlNI  GO  TO  20 

91. 

00 

NGES=NGES*1 

99. 

00 

IF  IN6ES.GT.MXGS1  GO  10  2C0 

9S. 

00 

lGEStNGES.NBlII 

96. 

00 

J6FSINGES.NB ISJ 

9T. 

00 

2C 

if  inu.eo.i.*no.j.ne.jmaxi  go  to  90 

98. 

00 

nchk=nchk«i 

99. 

00 

1CHk»nchk,nbi=I 

SO. 

00 

jchkinChk.nbi:j 

SI. 

00 

9  C 

CONTINUE 

S2« 

00 

NEVtNBl=NGCS 

SI. 

00 

IF  INB.EO.NBLK 1  60  TO  100 

S9. 

on 

NLNK=0 

ss. 

00 

00  80  K=1,NGES 

S6. 

00 

1=1CHK(K,N8> 

ST. 

on 

J  =  JCHK IK ,N8  1 

S8. 

00 

IF  INTVPEI l.J'll.EO.O)  GO  TO  80 

S9  • 

00 

NlNK=NLNK«1 

60. 

on 

L INK  1 NLNK ,NB 1=K 

61. 

00 

8C 

CONTINUE 

62. 

00 

NLINK INB 1=NLNK 

63  • 

00 

hxlnkihaxjihxlnk .NLNK9NGESI 

66  • 

00 

LllNB.ll=LltNBI*NGFS**2 

6$  • 

00 

IF  INB.NE.NBLKll  L 2 1 NB« I 1=L2 t NB l*NGES*NLNK 

66  • 

00 

LRlNVl=L2INBl*NGrS*NLN« 

67* 

00 

ICO 

JNIN=JHAX*1 

68. 

00 

C  IF  MXLNK.GT.NGES9NGES,  THERE  IS  SOME  WASTEO  STORAGE  THAT  IS  VERT 

69  . 

00 

C  OIFFJCULT  10  PROGRAM  AROOMO,  BECAUSE  OF  THE  DUAL  ROLT  OF  THE  LAST 

70. 

00 

C  BLOCK  OF  R1MV,  WHICH  IS  USEO  AS  A  SCRATCH  ARE*  FOR  RINV1  DATA. 

71  . 

00 

LRINV=L1INBLK 1 *H  A  X  0 1 MXLNK ,NGES**2l 

7?. 

00 

WR1TEI8.1S01  LRINV.LRINVI 

73. 

00 

ISO 

FOPHATtSSH  DIMENSIONS  RE0U1RE0  FOR  RINV  ANO  RINVl  ART  L01,  L02  =  , 

76. 

00 

A1S,2X,IS,1h./12uh  IF  EITHER  INPUT  VALUr  L01  OR  L02  IS  TOO  SHALL,  T 

75. 

oo 

BltE  COMPUTATION  IS  HALTEO  HrRr,  ANO  THE  USER  must  INCREASE  L01  ANO/ 

76  . 

00 

COR  L02  T0/22M  ThT  INOlCATEO  VALUfS.l 

77. 

00 

IF  HLOl. IT. LRINV1. OP. IL02.LT. LR1NV11)  STOP 

78. 

00 

RETURN 

79. 

00 

2  CO 

WRITE  16,202) 

80. 

00 

2  f? 

F0RMATI92M  OIMENSION  PARAMETER  MXGS  IS  TOO  SHALL  FOR  SPEC1F1C0  TOP 

81. 

00 

AOGBAPht  ANO  SUBREGION  SPECIFICATION.) 

82  . 

00 

STOP 

83. 

00 

End 
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VERSION  2 
(continued) 


SUBROUTINE  NtTINVIB.N.NI 
01NENSI0N  BIN,]|,B1I100),B2I100I 
NUN-1 

00  110  1=1, N1 
imm./Bii.n 
Bll, 11=1.0 
00  112  J=I.N 
112  6II,JI=BII,Jl*81tll 
121=1*1 

00  120  I1=IP1,N 
120  B1II1I=BII1,II 
00  12S  1 1  =  IP1 ,N 
125  Bill. 11=0. 

00  127  J=1,N 
127  621  Jl  =B  1 1 .  J I 

00  135  1 1  =  IP  1 *N 
DO  135  J=1.N 

135  Blll,JI=6IIl,JI-61IIll*B2IJI 
110  CONTINUE 

61111=1. /B IN, N1 
B IN.Niri. 

00  1*0  J=1.N 
1*0  BIN,J1=BIH,J1*61I11 

00  150  1=2, N 
00  155  12=1,1 
1*5  Bll 12  I =B 1 12, 1 1 
INI =1 -1 

00  156  12=1, INI 
1*6  6112,11=0. 

00  157  J=1 ,N 
1*7  B2IJI=BI1,JI 
INI = I —1 

00  160  12=1, INI 
00  160  J=1,N 

160  BII2,JI=BII2,J1-B 11121 *B2IJ1 
1*0  CONTINUE 


VERSION  3 


PR06RAH  JAYNORIINPUT,  OUTPUT,  T APE S : I NPUl , T A  PE  6 rOUTPUT  I 
C0HH0N/SEVP/AXI7,9I,Ay<7,9I,BBI7,9I,CXI7,9I,Q|7,9I,FI7,9I, 

A  RE  CUR  I  9,1 1 1  ,  R I N  V  I  S 68  I  ,  R I N  V 1 1  308  I  ,OUHOl  1»  ,21  ,OUHI  1191  • 

B  <1UH7<I<II  ,X|9, 11 1  ,H|9, 2  I  I.ERRtB.Ill.XXIB.m.ILFT  111*. 

C  IRGT'Iii,I6ES'I9,3I,JGES'I9,3I,ICHMI9,3i,JChm19,3I,NEVI3i, 

0  LINK|7,2I,NLINIH2I,NTYPEI9,III,IEI1I,LII31,L2I2I 
ORTA  H,N,HXGS,NBLK,L0I,L02/9,I1,19,3,S66,30B/ 

OATA  IE/9,7,10/ 

C  OIETRICH’S  H001FIE0  HAOALA  SOLVER,  VERSION  3. 

C  INIS  VERSION  IS  20ENT1CAL  TO  VERSION  2,  EXCEPT 

C  CY  HAS  BEEN  EL1H1NATE0,  < BY  01V101NG  THROUGH  THE  ELLIPTIC  EQUATION 

C  BY  THE  COEFFICIENT  OF  Xtl»l,J»2l,  THE  USER  CAN  USUALLY  GET  THIS  FORM 

C  IN  A  STRAIGHTFORWARD  HANNER,  IF  THIS  COEFFICIENT  VANISHES  NEAR  UPPER 

C  J-BOUNOAR IE S ,  A  H00IF1E0  FORK  OF  THIS  VERSION,  IN  WHICH  THE  RESIDUAL 

c  calculation  at  the  corresponding  error  points  is  hooifieo,  moulo  be 

C  USEFUL. I 

C  EQUATION  TO  BE  SOLVEO  IS 

C  NTYPE  I  1M,  J*I  I  •  lAXII.JI  •  XII.J'll  «  CXII.Jl  p  XII«2,J«I> 

C  «  AYll.JI  9  XII'I.JI  ♦  X  < I « 1 , J«2 I  I 

C  ♦  12  -  NTYPE  II«I,J«1 11  p  BBU.Jl  P  XII«1,J*1I 

C  :  2  P  II  •  NTYPE 1 1«I, J*1 1 1  P  BBI I ,JI  P  XII*I,J*II 

C  «  NTYPEII*I,J*II  •  F  < I , Jl ,  I  1 1 :0,I, . .. ,H-1 1 , I JtO.l , . .. ,N-1 1 1 

C  WHERE  TERRS  OUTSIOE  OIHENSION  BOUNDS  ARE  INTERPRETED  AS  2ER0*S. 

C  THUS,  IF  NTYPEl I,J|:0,  XII, Jl  IS  LEFT  UNCHAN6E0, 

HltH-I 

NItN-I 

H2tH-2 

N2:N-2 

NBLKltNBLK-I 

00  BO  J=1,N 

oo  an  i=i, h 

sr  NTYPEl I, Jl=0. 

00  90  J:2,N1 
00  9(1  1:2, HI 
9T  NTYPE II,Ji:l 

00  2000  JCt2.Nl 
00  2000  ICt2,Hl 
NTYPE IIC, JCItO 
NTYPEl IC-I.JCISO 
NTYPEl IC«I, JCItO 
NTYPEIIC.JC-lltO 
NTYPElIC, JC«I|:0 

WRITE  I  6,92  I  1 1 NTYPE 11 , Jl ,1 1 1 ,H I , Jtl ,N I 
92  F0»HATII0Xf9III 

CALL  T0P0LIRECUR,1LFT,IRGT,IGES,JGES,ICHN,JCHN,NEV,LINK,NLINK, 

A  NTYPE,IE,H,N,H1,N1 ,H2,NBLK,NBLKI,HXGS,LI,L2,L0I,L02I 
00  100  Jtl,N2 
00  100  1 1 1  ,H  2 
AX  I I,JltI.*.OIP| I-9IPP2 
CXII,JI:I.-.0IP|I-9|PP2 
AY  1 1 ,Jlt 1 . y.OIPI J-5IPP2 

BBII,Ji:i*AX<I,JI-AYI2,JI*CXtI,JI*.0lPIJ*SIPP2*I.II/ll.*.0lPIJ*5l 

APP2I 

100  CONTINUE 

OENOHtl .2SPHPNIPP2 
00  170  Jt 1 ,N 
00  170  II1.H 


VERSION  3 
(cont inued) 


XI 1,  Ji:i  1-1.  I9I1-M  tOj-l.  1*1  J-NI/OENOH 
XXIl,JI=Xll,J|9|l.-NTYPE)l,J)l 
170  CONTINUE 
FSUN=0. 

00  171  J=2,N1 
l=J-l 

00  171  1=2, N1 
K  =  l-1 

F)K,LI=AXIK,L|9X)l,JI«AYIK,Ll9X)l,LI*CXIK,LI*Xll«l,J>* 

A  XI1,J*1I«BBIK,l1*X|1,JI 

i7i  fsum=fsun»absifix,lii 

FSUN=N?*N2/FSUH 
T0=SEC0N0l TIN ) 

CALL  JAVIAX.AV^B.CX.RINV.RINVI.H, RECUR, 1LFT,IR6T,1GES,J6ES, 
A1CHX,JCHX,NEV,L1NK,NL1NX,NTYPE,1E,N,N,M2,N2,NBLK,NBL*1,HX6S,L1,L2I 
T1MSUN=0. 

00  1B1  J=2,N1 
00  181  1=2. HI 
011-1 ,J-1 I  =F  I 1-ltJ-ll 
1(1  Xll,JI=0. 

ALINE=0. 

7 0  =  SEC0N0l TIN  I 

CALL  NOR  1AX ,AY , BB , CX , R 1NV , R IN  V 1 ,0  UM0,0UH1 , OU*  2 ,0,H, XX, RECUR, 

A  1LFT,1RGT,1GES,JGES,1CHK,JCHK,NEV,L1NK,NL1NK,1E,M,N,M2,N2,NBLX« 

B  N6L*1,MXGS,L1«L21 
00  990  J=2.N1 
00  990  1=2, N1 

990  XI1,JI=XI1,JI«XXI1,JI 
00  991  J=1,N 

00  991  1=1, N 

991  XXI1,JI=0. 

00  992  J=2,N1 
L=J-1 

00  992  1=2.N1 
*=1-1 

992  OI«,LI=FIK,LI-AYIK,L>*Xll,J-ll-AX|K,L>9Xll-l,J>-BB>K,L>*Xll,JI- 
A  CXIK,L>»Xll’l,JI-Xll,J,l  ) 

1  COO  CONTINUE 

9  C5  F0RNAII1X, 1P1GEB.11 
ERSUN=0. 

00  790  1=2, N1 
*=1-1 

00  7*0  J=2,N1 
L=J-1 

errii,ji=fik,l)-ayix,l>*xii, j-1 i-axi*,li*xii-i,ji-bbi*,l>* 

A  Xll,JI-CXIK,LI9Xtl.]vJ|-Xll,J*l) 

ERRI1,JI=ERRI1,J1*NTYPEI1,J1 
ERSUM=ERSUN**BSIERRll,JII 
790  CRRI1,JI=ERRI1,UI9FSUN 

ERSUH=ERSUH*FSUH/IN29N2 1 
WRITE  IB, 9091 

9C9  F0RNATI/90X.2OHN0RNAL12ED  RESIOUAlSI 
00  792  1=2, N1 

VRlTEIt,905M£RRIl,JI,J  =  2,Nl> 

792  CONTINUE 

WRITE  16,906  I  ERSUM 

9 CG  F0RNATI/J7H  NEAN  ABSOLUTE  N0RNAL17E0  RES10UAL  =  .1PE9.2I 
IF  IERSUN.GT.1.E-9I  STOP 
IF  I1C.NE.2I  NTYPE11C-1.JC)=1 
IF  IlC. NE. Mil  NTYPEI1C«1,JCI=1 
IF  IJC.NE.2)  NTYPEIlC,JC-l)=l 
IF  IJC.NE.N1I  NTYPEI1C,JC*11=1 

2  COO  NTYPE 1 1C,JC)=1 


VERSION  3 
(continued) 


!• 

01 

subroutine  j«m>,i»,siiCx ,rinv, rinvi.h, recur. ilft. 

2« 

00 

1 1RGT «1G£S • JGE S ,  1 CHM ,  JC  HK ,  NE  V .LINK .NL Ink »NTTPE  v 1E»hO,NO,M2, 

3. 

00 

*  NBLK,NBLKl,HXGS,Ll,L2l 

«• 

00 

01  PENSION  XX (H2.N2 1 ,IT)h2,N2),8B(H2,N2I ,C» IH2.N2I , 

5. 

00 

*  R1NVI1I.R1NVK1  >, REC UP l"0, NOl.HIHO.NOl, 

6  • 

00 

6  ILFT(NO>,lRGT(NOl ,L1(NBLKI,l2(NBLK1I,16ES(HXG5,NBLK>, 

?. 

00 

C  J6ESfHXGS,N8LKI,ICHK(HXGStNBLK  1,  JCHX IMXGS ,NBLK 1  .NEVtNBLXI 

«. 

00 

0  LINM(H2,NBLX1 I.NL1NK INBLK1 1 tNTTPE 1H0.N0 1 , IE INBLR I 

9. 

00 

JL=1 

10. 

00 

L  OIL  1 1 NBLK 1 

n. 

00 

NB=0 

12. 

00 

iro 

NBINB* 1 

11. 

00 

NGE  SINE  V INB 1 

m. 

00 

NGESIiNEVINB-1 1 

is. 

00 

LAiLl(NB) 

lto  • 

00 

LBIL2INB-1 1 

IT. 

00 

LC-L2I NB 1 

IS. 

00 

JH-IElNB t 

19. 

00 

JHPIJHM 

re. 

00 

JHMUH-2 

21 . 

00 

00  253  NGi 1 . NGES 

22. 

00 

IGIIGESING.NB  1 

23. 

00 

JG- JGE  S (NG . NB ) 

29. 

00 

JF-JG-I 

rs. 

00 

00  210  J- JL « JHP 

2b. 

00 

00  210  1=1, HO 

27. 

00 

210 

Mil , Jl=0. 

28. 

00 

HUG,  JG  1  =  1. 

29. 

00 

IF  IJF.GT.JHMI  GO  70  228 

10. 

00 

IF  INTTPEIlG.JFj.EO.Ol  GO  TO  220 

31. 

00 

NOPiNLINX INB-1 1 

12. 

00 

00  219  Nil, NOP 

33. 

00 

HILIN* IN.NB-ll 

39. 

00 

IF  (ICHxiH.NP-lI.EC.lG.ANO.JCMHIM.NB-ll.EO.JFI  GO  TO  215 

IS. 

00 

219 

CONTINUE 

3b. 

00 

215 

00  218  Nil  * N OP 

37. 

CO 

LILINMN.NB-ll 

38. 

00 

2  IS 

HtICHKtL,NP-ll,JF|IRlNYl)LB«)N-l >*N6ES1«MI 

39. 

00 

2  20 

00  225  JIJF.JHM 

90. 

00 

ILHLFTIJ.2I-I 

91. 

00 

1RH*GT|  J.21-1 

92. 

CO 

2  22 

00  225  HIl.lR 

93. 

00 

2  25 

H(l«l,  J«2li-RECURI1«1,J.2>*I1XII,J>*H)I,J»I|.AVII,JI*HII*1 

99. 

00 

A  BBII, JI*H(I41,J«1 t«cx)I,JI*H|l«2,J*I) t 

95. 

00 

2  28 

oo  230  nch.nges 

9b. 

00 

1 IIC  HK INC.NS l-I 

97. 

00 

JiJCHKINC.NBI-l 

96. 

00 

2  30 

RINVILX.(NC-1|*NGE5*NGIIXXI1  ,JI*H)I,J«I).ATI  1,JI*HI  LI,  J|.| 

99. 

00 

a  hu«i,j»i t.cx(i,j)«Mii.r, j«i> 

SO. 

00 

IF  INB.EO.NBLKI  GO  TO  250 

SI. 

00 

NOPINLlNKiNB) 

S2. 

00 

JHEINBI 

S3. 

00 

00  29C  Nil, NOP 

$9. 

00 

MILINX  1N.NB 1 

ss. 

00 

2  90 

RINVILOON-l  l*NGE5«NG>IHI  ICHK(H,NB>,  ji 

Sb. 

00 

2  30 

HI  16, JG!=0. 

57. 

00 

CALL  HI  TINVIRlNVlLX.il, NGE5, NGE S» 

Sb. 

00 

IF  INB. EO. NBLK I  RETURN 

59. 

00 

00  260  IH. NGES 

SO. 

00 

00  260  Jil.NOP 

bl. 

00 

R1NVULC«IJ-I)9NGE5.IU0. 

82. 

00 

00  260  KI1.NGES 

S3. 

00 

2(0 

•1NV1ILC* )J>1I*NGE5«I 1IR1NV1 ILC.lJ-1 )*NGES«1 i - 

89. 

00 

A  RINTILA.IK-I  I9NGE5.I I9R1NVIL0. I J-l I9NGES.R 1 

SS. 

00 

C  R 1 NV 1(1 ,  J ,  NB 1  IS  THE  ,,«l*OST  HOMOGENEOUS*  *  S0LU1I0N  »1  the  j 

bb. 

00 

C  RTS10UAL  P05IT1  ON  FORCEO  S3  •  RE510UAL  VSIUE  OF  1  »T  THE  I-lH 

87. 

00 

C  Rf S 1 OU A  L  POSITION.  H0H0GENE0U5  B.C.’S  IRE  AS5UHE0  E TERTNHERE i 

SS. 

00 

c  iNCLuoiN6  the  top  of  the  present  susregion  ns. 

89. 

CO 

JLIJH 

70. 

00 

CO  10  100 

71. 

00 

END 

29 


VERSION  3 
(continued) 


1. 

01 

SUB40UT1NE  NOPIAX,AX,B6,CX,PlN9.41NVl.OUHO.OUH1.0UH'Z,F,H,X, 

2. 

00 

IPECUR, 1LFT.1RGT, 1GES. J6E S , 1 CHK , JCMK ,NE V .LINK .NL1NK ,1E .HO.N0.N2, 

3. 

00 

*  N2.N8LK.NeLKl.8XGS.Ll.L2t 

4. 

00 

0 I NENSI ON  AXIH2,N2l,AVIH2,N2l,8BIH2,N2l,CXIN2,N2l, 

S. 

00 

4  9INVIll,RlN91lll,9tCURIN0,Nnt,HIHD,N0l, 

4. 

00 

B  ILFTINOt .IRGTlNO 1 ,L1 INBLK I ,L2 INBLK 1 1 ,1GES IHXGS ,NBLK I , 

7. 

00 

C  J6E S 1 MXGS iNBLK 1, 1 CNK 1 HXG S.NB LK 1, JCHK (HXGS.NBLKI.NE 9 INBLK 1, 

4. 

00 

0  LINK l«2. NBLKII.NL INK INBLK 11 .IE INBLK 1 

9. 

00 

01HENS10N  0UN0IH2,NBLKll,0UHl|HXGSt ,0UN2 IHXGS 1 ,F IM2.N2 1 ,X I  HO, NO  1 

10. 

00 

jsn 

11. 

00 

00  ISO  N8:l.NBLK 

12. 

00 

LCXL2INBI 

13. 

00 

JFllEINBt-2 

14. 

00 

00  IPS  JXJS.JF 

IS. 

00 

lLxitFTI J«7l-1 

14. 

00 

IArI»6T|J.2l-I 

17. 

00 

00  IDS  1 -1L , 19 

14. 

00 

ICS 

XII«l,J«2tXPECUPtl«I,J*2l4IFIl,JI*AXII,JI4XII,J.lfAVIl,J|4 

19. 

00 

A  XI1M,J>-BBII,JI4X|I«I,J*1I-CXII,J>4XI]*2,J«]II<I1.- 

20. 

00 

B  9ECU9I1*1,J«2II*XI1»1,J«2> 

21. 

00 

IF  INB.EO.NBLKt  60  70  ISO 

22. 

00 

NGES:NC9INBI 

23. 

00 

00  115  NX  1 , NGE S 

24. 

00 

IxiCHKIN.NBI-l 

2S. 

00 

JXJCHKIN.NBI-l 

24. 

00 

1  IS 

0UH1 INIXF 1 I , JI-AXI 1 , J|*XI I ,J«1 l-AT 1 I , JI4XI I«1 , JI-BB I  I, Jl* 

27. 

00 

4  *ll»l,J»lt-CXII,JtAXI1.2,J»l>-XIIM,J*2l 

24. 

00 

NOPxNL INK  1 NB  t 

29. 

00 

JX1EINBI 

30. 

00 

00  120  NX], NOP 

31. 

00 

0UH2INIX0. 

32. 

00 

00  118  NX1.N6ES 

33. 

00 

114 

0UH2INH0U«2  INI  *0081 IHI491N91ILC*  IN-1  I^NGES^I 

34. 

00 

nxlinkin.nbi 

3S. 

00 

I X I CHK 1 H ,NB 1 

34. 

00 

OUHOIN.NB txxil.J) 

37. 

00 

1  20 

XI1,J|XXII,JI-0UH2INI 

34. 

00 

ISO 

JSXIFINBI 

39. 

00 

00  300  NB1=1,NBLK 

40. 

00 

NBXNBLK-NB1«1 

41. 

00 

JStI 

42. 

00 

IF  INB.NE.II  JSxlEINB-II 

43. 

00 

JFxiE 1 NB 1 -2 

44. 

00 

LAXLllNSt 

45. 

00 

IBXL7INB-1 1 

44. 

00 

NGE5XNE4 INB ) 

47. 

00 

IF  INB.EO.NBLKt  60  70  201 

48. 

00 

JxIEINBI 

49. 

00 

NOPXNLINKINBI 

SO. 

00 

00  200  NX], NOP 

SI. 

00 

H XL  INKIN, NB 1 

52. 

00 

IHChkih.nb  i 

S3. 

00 

2  CO 

Xll.JHOUHriN.NBt 

54. 

00 

2  Cl 

CONTINUE 

■s. 

00 

NX1EINBI 

54. 

00 

00  202  JXJS.N 

S7  • 

00 

00  202  I XI .NO 

30 


VERSION  3 
(continued) 


2C2  HII.JCD. 

00  210  N=1.N6£S 
I=ICHKIN.NP)-1 
J=JCHK IN.NBI-1 

2  10  0uHHNlrriI,JI-AX(I,Jt*XtI,J*l«-ATtI,JI*Xtl»I,JI-66ll,JI* 

«  X(l«l.JM(-CX(ltJ(*X(l«2,J*l(-Xll*ltJ*2l 
00  220  N-l f M6E S 
DUN2INU0. 

DO  218  MM.NSES 

218  0UH2IN>I0UH2INI«0UHHH«*RlN¥tLA»IN-ll*N6ES»HI 
I=1GES(N,NBI 
JTJGfSIN.NBI 
HI  2  . JI=DUH2|N) 

220  Xll.JlxxiI.jl ♦ DUH2 (Ml 
IF  INB.E4.il  GO  TO  2SD 
DO  222  N=1,MXGS 
222  OUNHNIXD. 

nop=nlinmnb*ii 

J:JS 

00  2 TO  N=1,N0P 
NIIINKIN.N8-1I 
I:ICHKItt.N6-ll 
210  0UN1<HUH(1,J«1I 

C  WHEN  POUNOARt  IS  IRREGULAR,  THERE  IS  SDHE  UASTEO  CALCULATION  HERE  IN 
C  ORDER  TO  A VL 1 0  LOGICAL  DECISIONS  OR  EXTRA  STORAGE  (NGES  COULD  BE 
C  REPLACED  Bt  NL1NKINBII 
NGES-NCVINB-1 I 
00  2«0  N-l .NOP 
0UN2(Nl:0. 

00  238  M - 1 VNGE  S 

2  38  0UH2INU0U"2<NI«0UH1(M<*R1NV1(LB*  IN-1I*NGES*H) 

H:lINK(N,NB*1 I 

2  AO  HI ICHK IH.NB-l I , JCMK (N.NB- 1  I  I =DUN2 IN  I 
2*0  DO  jro  JXJS.JF 

ILXILrTI J.21-1 
1R-1RGT I J. 2 | -1 
00  300  KIL.IR 

0UHl(ni-RECURU»».J»21*l»X(l,J»*HII,J*ll«ATII,JI»H(I*l,JI» 

A  BBII,J(*HI1«1 ,J«ll»CXIltJ|AM|I.2,J.l(l 
H(I«1  tj.?i:DUHli  !(<■  1.-RECURI  l*l.J*2M*Hll*l,j*2l 

3  CO  XI  1*1,  J*2CIU<|,J.2l*0u»lUI 


VERSION  3 
(continued) 


1. 

01 

SUBROUTINE  T0P0LIRECUR,ILFT,IRGT,IGEStJ6EStICHK,JCHK,NE¥tLINK, 

7. 

00 

1NLINK.NTYPE.IE  ,M,N,Hl,Nl,l«?,NBlll,NBlKl,HX&S,Ll,l?,l01,L0?l 

1. 

00 

OIMENSTON  RECURIP.NI,  IlF T I N I , I RGT I N I , IGE S I HXG S .NBLK I , 

4. 

00 

*  UGE  SIHXGSfNBLX 1 t IC  MX IHXGS.NBL* I «  JCHK  IHXGSiNBLK I . N  E  V 1 NBLK I , 

s. 

00 

6LINKIH?vNBLKl).NLINKINSLKltvNTYPEtH,Nt,IEtNBLKtvLltN6LKt,L?tN6LKl) 

i. 

00 

00  10  J=?,N1 

T. 

00 

1=0 

a. 

00 

7 

1=1*1 

9. 

00 

IF  tNTYPEtl,Jt*NT¥PCtI,J-ll.EO.O.*NO.I.LT.Hl>  GO  10  7 

10. 

CO 

Ilftij»=i 

11. 

00 

I=h 

12. 

00 

4 

1  =  1-1 

13. 

00 

IF  INTYPEtI,J>RNTVPEfItJ-ll.E0.0.*N0.I.6T.l>  go  to  4 

14. 

00 

10 

IRGTt Jt=I 

IS. 

00 

IRGTI 11=7 

1  B  • 

00 

ILFTI1C? 

17. 

00 

ILFTINC? 

18. 

00 

1 RGT IN l=? 

19. 

00 

00  15  J=?,N1 

70. 

00 

IF  IILFTI JI.LE. IRGTI Jll  60  TO  15 

71. 

00 

IRGTI Jl=7 

77. 

00 

ILFTIJI=2 

73. 

00 

15 

CONTINUE 

74. 

00 

JMIN=? 

?S. 

00 

HXLNKZO 

?6  • 

00 

11911=0 

77. 

00 

L7I 1 1=0 

78. 

00 

00  100  NB=1,NBLK 

79. 

CO 

N6ES=0 

30. 

CO 

N  CHICO 

31. 

00 

JM»X=IE INBI 

37. 

00 

00  40  I =7 , M 1 

33. 

00 

nh=ntypeii,jhin-i» 

34. 

00 

NU=NTyPE9I.J"INI 

3S  . 

00 

00  40  J:JH1n(JN6X 

36. 

00 

nl=np 

37. 

00 

NH=NU 

38. 

00 

RECUR  1 1 . J l=NH 

39. 

00 

IF  INL.EO.OI  RECUR  1 1 .J 1=0. 

40. 

00 

NU=NTYPEII,J*1I 

41. 

00 

IF  INH.EO.OI  60  TO  40 

42. 

CO 

IF  INL.E0.1.6NO.J.NE.UHINI  GO  TO  70 

43. 

00 

NGE  S=NGE  5* 1 

44. 

00 

IF  INGE5 .GT • HXG5 I  GO  TO  700 

45. 

00 

IGESINGES.NBCI 

46. 

CO 

JGE5INGE5.N6CJ 

47. 

00 

70 

IF  INU.EO. I • XNO.J.NE • UHX X 1  GO  TO  40 

48. 

00 

NCHICNCHK.l 

49. 

00 

ichk INCHK.NB 1=1 

50. 

00 

JCHK INCHK.NB  »=j 

SI. 

00 

40 

CONTINUE 

57. 

00 

NEVINB  »=NGE5 

S3. 

00 

IF  9N8.EO.NBLK)  GO  TO  100 

54  . 

00 

NINICO 

SS. 

00 

00  BO  4=1 tNGES 

56. 

00 

I  =1 C MX  IK  .NB  1 

57. 

00 

J:JCHKIK,N6> 

58. 

00 

IF  INTYPEII.J«1I.EO.OI  GO  TO  80 

59. 

00 

NLNK=NLNK«1 

60. 

00 

IINKINLNK.NBCK 

61. 

00 

80 

CONTINUE 

67. 

00 

NLINKINBCNLNK 

63. 

00 

HXLNK=««XOIHXl.NK  .NLNK4NGE5) 

M  • 

00 

11INEM  t=LI  INBI'NGES**? 

65. 

00 

IF  9N6.NE.N6LK I  1  LlINB.I l=L7INB (»NGE5*NLNK 

66. 

00 

LRIN¥1=L?INB>*NGES*NLNK 

67. 

00 

100 

JHlN=JH6X*I 

68. 

00 

C  IF 

NXLNK.GT.NGE54NGE5.  THERE  IS  50HE  H85IE0  5T0R8GE  TH»T  IS  XERY 

69. 

00 

c  oifficult  to  rrogra*  arouno*  because  of  the  OUAL  ROLE  of  the  last 

70. 

00 

C  BLOCK  OF  ftlNVt  WHICH  IS  USEO  AS  A  SCRATCH  AREA  FOR  RINVl  0 ATA  • 

71. 

00 

LR  IN  V  =L  J  INBLX)«M»xOl><XLNK  .NOES**?! 

77. 

00 

WRITE  16.1501  IPINX.IPINVI 

73. 

00 

150 

F0PH6TIS5M  OIHENSIONS  REOUIREO  FOR  RlNV  *N0  RINYI  »RE  L01,  LD2  =  . 

74. 

00 

6IS.7X.1S.1H./174H  IF  either  INPUT  »»IUE  LOI  OR  LO?  15  TOO  SH8LL •  T 

75. 

00 

BHE  COMPUTATION  IS  H4LTE0  HERE,  »N0  the  USER  "U5T  INCRE85E  LOI  »N0/ 

76. 

00 

COR  LO?  TO/72H  THE  IN0IC8TE0  ¥»LUCS.» 

77. 

00 

IF  99L01.LT.LRIN¥».0».IL07.LT.IRIN¥1>>  stop 

78. 

00 

return 

79. 

0” 

7  CO 

WR 1 TE  9  6, 7071 

80. 

00 

7  C? 

F0RH6TI97H  01HEN510N  P6RXHETER  HXG5  IS  TOO  S"»LL  FOR  SPEC1FIE0  TOP 

81. 

00 

60GR6PHY  ino  subregion  specif icmion.i 

87. 

00 

stop 

83. 

00 

ENO 
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VERSION  3 
(continued) 


SUBROUTINE  MlTlNVlBtN.NI 
DIMENSION  B  I  N  ,  1 I .B1 1 1  00  I ,  B  2 1 1001 
N1IN-1 

00  110  1*1 .Ml 
Bll 11 =1 ./Bl 1 .1  I 
Bll, 11=1.0 
00  112  J=1 «M 
11?  B<1,JI=B<1,JI*B111I 

1P1=1*1 

00  120  11=1P1,M 
120  B 1  <  1 1 1  — B  « 1 1 1 1 1 
00  12S  11=1P1,N 
1  25  Bill. ICO. 

00  127  J=1,N 
127  B2IJ)=B«1,JI 

00  135  11=1P1,M 
00  135  J- 1 « M 

135  B<ll,JI=BIIl,JI*Blllll*B?IJt 
110  CONTINUE 

B1<1I=1./BIM,MI 
8  IH,N)=1. 

00  110  J-l t M 
110  B|M,JCBIM,JI*B1I1I 
00  150  1=2. N 
00  155  12=1,1 

1 55  81<12I=BI12,1I 
INI =1*1 

00  155  12=1, INI 

156  8112,11  =  0. 

00  157  J=1,H 

157  B?IJI=SI1,JI 

lMlil-1 

00  160  17=1, INI 
00  160  J=1,H 

160  B< 12,JI  =  B< 12, Jl-Bll  12i*B?IJl 
150  CONTINUE 


VERSION  4 


»  • 

r. 


i. 

00 

PROGRAM  MAIN) INPUT, OUTPUT , TAPE  5  =  1 NPUT ,T  APE4-0UTPUT ) 

2. 

00 

C0MH0N/SEvP/AX(7,4) ,AYI7,4t ,6617,41 ,CX(7,4(,0(7,4I,F(7,4(, 

3. 

00 

A  RlNVlSAA>,RlNVI(304l,OUHO(14,2(,0UM]l]4t, 

4. 

00 

4  OUM2tl«t,X(4,lII,Ht*,]lt ,ERR|4,11I ,XXI4,Ilt,ILFTtlII, 

s. 

00 

C  IRGTI III.IGES 114,31, JGESt 14,31 .ICHKI 14,31, JCHK I 14,3t,NEVI3 I, 

4. 

00 

0  LlNKI7,2l,NLINKI2l,NTYPE t4,lIt,IEI3t,LII3t,L2t2> 

7. 

00 

OAT A  M,N,HXGS,N6LK,LQI,L02/4,I1,14,3»546,306/ 

(. 

00 

OATA  IE/4, 7, 10/ 

(. 

00 

c 

oietrich’s  hooifieo  maoala  solver ,  version  4. 

id. 

00 

c 

THIS  VERSION  is  IOENIICAl  10  VERSION  3,  EXCEPT 

ii. 

00 

c 

NO  ISLAN05  ARE  ALLOUEO  ANO  IRREGULAR  4OUN0ARIE5  MUST  BE  SUCH  THAT  ALL 

12. 

00 

c 

N0N2ER0  NTYPE  VALUES  IN  A  GIVEN  RON  MUST  BE  CONTIGUOUS,  THEREBY 

13. 

00 

c 

requiring  that  the  right  and  lefi  boundary  curves  be  single  valued 

14. 

00 

c 

FUNCTIONS  OF  Y.  THIS  ELIMINATES  THE  NEEO  FOR  ARRAY  RECUR  ANO 

1$. 

00 

c 

associateo  computations. 

% 

II 

14. 

00 

c 

EQUATION  TO  BE  SOLVED  IS 

17. 

00 

c 

NTYPE(1*1,J«I I  •  (AXII.J)  •  Xtl,J*It  «  CXII.JI  •  X ( I . 2 , J.l 1 

14. 

00 

c 

«  A  Y  1 1  ,  J )  •  XtlM.Jt  «  X|I«1,J»2H 

It. 

00 

c 

♦  (2  -  NTYPE 1 1«1 , J«1 1 1  •  46) I • J(  •  X)I.],J.]| 

20. 

00 

c 

r  2  •  11  -  NTYPE II» I, J«1 1)  •  46(1 , J 1  •  XII.l.J.l) 

21. 

00 

c 

«  NTYPE(I«1,J«1(  •  FII.JI,  t(I=0,I,...»H-It,IJ=0,l,...,N'l)t 

22. 

00 

c 

HHERE  IERMS  OUTSIDE  DIMENSION  BOUNOS  ARE  INTERPRETED  AS  2ER0*S. 

23. 

00 

c 

THUS,  IF  NTYPEII,  JCO.  XII.JI  IS  LEFT  UNCHANGED. 

24. 

00 

MI=M-1 

'  • 

25. 

00 

N1=N-1 

24. 

CO 

M2=M-2 

-  H 

27. 

00 

N2-N-2 

24. 

00 

N6LKl=NBLK-t 

2«. 

00 

00  40  J=1 ,N 

30. 

00 

00  40  1=1, M 

31. 

00 

4  C  NTYPEII, Jl=0. 

32. 

00 

00  40  J=2.N1 

33. 

00 

00  40  1=2, MI 

34. 

00 

4  C  NTTPm.JCl 

35. 

01 

NTYPE 12, 31 .NTYPE 12,4 », NTYPE 13,4 1, NTYPE 12,5 1, NTYPE 13,51 .NTYPE 12,4) 

34. 

00 

A  =0 

• 

37. 

00 

WRITE  14, 421  IINTYPEIt,JI,t  =  I,HI,J=l,NI 

■t 

.  fc 

34. 

00 

42  FORUM  liox.411 1 

34. 

02 

CALL  TOPOL  1  ILF T, IRGT , I6ES , JGE S , I CHK , JCHK ,NE V . LINK, NLINK, 

. , 

40. 

00 

A  NTYPE, IE, M,N,MI»N1,M2, NBLM , NBLK 1 , M XGS , L 1 ,L2 , LOI ,L02I 

41. 

00 

00  ICO  J=I,N2 

r. 

42. 

00 

00  100  1=1, M2 

Lgh 

43. 

00 

AXII.J 1=1. ..01*11-41 **2 

\9 

44. 

00 

CXI1,JI=1.-.0I*)I-4I**2 

M 

4$. 

00 

AyI1,J(=1.«.01*IJ-SI**2 

|  „*  • 

44. 

00 

4BI1 ,JI=I-AXII .JI-AYII , J 1 -CX 1 I ,J)..0I*IJ-SI**2-I. ll/ll.-. 01*1 J- 51 

*4 

L 

47. 

00 

A**2| 

!*• . 

44. 

CO 

100  CONI INUE 

44. 

00 

OENO“=l.2S*M*NI**2 

50. 

00 

00  170  J=1,N 

J 

51. 

00 

00  170  I=I,M 

52. 

00 

XII, Jinl-l. 1*11 -M  1*1  J-I,  )  *1  J-N  )  /  OtNOM 

* 

• 

53. 

00 

xxii,ji=xii,ji*ii.-ntvpeii*jm 

54. 

00 

170  CONTINUE 

N— 

t- 

55. 

00 

F  SUM  =  0 , 

54. 

00 

00  171  J=2,N1 

57. 

00 

L=J-1 

V 

f: 

= 

Lj 

• 

• 

' 

- 

> 
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VERSION  4 
(continued) 


sa* 

00 

00  171  1:2, Ml 

59. 

00 

«=I-I 

to* 

00 

FtK,LI=AXIK,L|9XI!,Ji«AVtl',LI*Xll,LI*CXIK,LI*XII*l,JI* 

tl* 

X  0 

*  III,J*1I'BBIK,LI*X|I,JI 

42. 

00 

I  71 

F$un=F$um«abSIFik,l»I 

€3. 

00 

FSUM=N2*N2/FSUM 

tb* 

00 

7  D-SFCONO 1 7 1 M  4 

65* 

02 

CALL  jAVIAX,Ar,BB,CX,RINV,PlNVl,M,lLFI,IR6T,I6ES,JGES, 

tb* 

00 

A ICMK fjCMX.NEV.LlNK ,NL 1 NK , N7 7 R E . IE ,M ,N , M2 ,N2, NBLK , NBLK 1 , NIBS. LI, L2 1 

00 

7 Im$uh:o. 

tB* 

CO 

00  181  J=2,NI 

49* 

00 

00  181  1:2. MI 

70* 

00 

011-1,  j-n=Fu-i,j-n 

7i. 

00 

181 

iii.jco. 

72. 

00 

ALIME-O. 

73* 

DO 

TO=StCONOlTlM| 

79  • 

02 

CALL  NOR IAX,Ay,86,CX,R!NV,PlNVl,DUN0,0UMl, OUM  2 , 0 , M , X  X , 

75. 

00 

A  ILFT.IRGI.IGES, JGE S.ICmk.jCMk , NEV.LINK.nl INK, IE ,M,N,M2,N2,NBLKf 

76. 

00 

B  NBLK1.MXGS.LI.L2I 

77. 

00 

00  990  J=2,NI 

7a  • 

00 

00  990  1^2, Ml 

79. 

00 

9  90 

xiI,ji=xii,jI«xxII,ji 

60. 

00 

00  991  j:  I  ,N 

61. 

00 

00  991  HI, M 

e:. 

00 

991 

XXI I,JI=0. 

as. 

00 

00  992  J=2,NI 

PH. 

00 

irj-l 

85. 

00 

00  992  1=2, MI 

at . 

00 

K  =  I-1 

.7. 

00 

9  92 

oik,li=fik,li-av«k,li*xii,j-ii-axik,li*xii-i,ji-bbik,l>*xii,ji> 

8k. 

00 

A  CXlK,LI9XII«I,Jt-XII,J»II 

17. 

00 

ItOQ 

CONTINUE 

*3. 

00 

9  CS 

FOPMAIIIX, IRISES. II 

*1. 

00 

ERSUMZO. 

*2. 

00 

00  7.0  I=2.MI 

93. 

00 

k  =  I-1 

... 

00 

00  7.0  J=2,M1 

»S. 

00 

t=J-l 

*6. 

00 

ERRlI,Jl=F  IK ,L I-ATIK.L I»X|I , J-l  >-AX|K,lt*X  1 1 -1  ,  J  > -BB  1 K  ,  L  >  • 

97. 

00 

A  XII,JI-CXIK,LI*XtI«I,JI-XII,J«II 

98. 

00 

ERRII,  JHERRIl  .Jl.NVVPE  II,  Jl 

99. 

00 

ERSUMlERSUM.ABSlERPlI.Jll 

ICO. 

00 

7.0 

£RBII,JI=E9RII,JI9FSUM 

101. 

00 

ERSUM;ERSUm»fsuM2IM28N2I 

102. 

00 

WRI7E 18 ,90. ) 

103. 

00 

9  C. 

FORMA] | /.O*. 20MNORMALI2E0  RESIDUALS! 

10.. 

00 

00  7.2  1=2, MI 

10$. 

00 

WPIIEI8.90SIIERRII,JI,J=2,N1I 

108. 

00 

7.2 

CONTINUE 

107. 

00 

WRITE  18,9061  EPSUM 

108. 

00 

9C8 

FOPMA7I/37M  ML  AN  ABSOLUTE  NORMALIZED  RESIDUAL  =  .IRE9.2I 

109. 

00 

IF  IFRSUM.G7.1.E-.7  SIOR 

ItO. 

00 

cno 

35 


VERSION  A 
(continued) 


1. 

01 

SUBROUTINE  JAYIAX.AY.BB.CX.RINV.RINVI.M.IlFT, 

2. 

00 

llRGItIGESfJGES.ICHK.JCHK ,NE  V .LINK ,NLlNK,NTYPE,IE ,M0 ,ND ,M2 , N2 , 

3. 

00 

A  NBLK,NBLK1,mxGS.L1.L2> 

4. 

00 

DIMENSION  AX(M2fN21,AY IM2 ,N2 » , BE 1 M2 ,N2 1 , CX IM2 ,N2 »  t 

s. 

00 

A  RlNVtll.RlNVltll.HIMO.NQI, 

6. 

00 

B  ILFTlNDl  f  IRGTINOI  ,L1  ( NBLK 1 ,L2 INBLK I 1 . IGES (HXGS ,NBLK  » , 

I. 

00 

C  JGESIMXGS.N6LK1 .ICHK (MXGS.NBLK 1 .JCHKIMXGS.NBLK I.NEVINBLKI, 

e. 

00 

0  LINMN2,NBLK1>  fNLlNK  INBLK  1  )  ,NT  YP  E  1  MO  ,N0  )  ,  IE  INBLK  1 

9. 

00 

JL-1 

ID. 

00 

LDxLlINBLKI 

11. 

00 

NBrO 

12. 

00 

ICO 

NBXNB«1 

13. 

00 

NGESXNEVINBl 

14. 

00 

NGES1XNEVINB-1 1 

IS. 

00 

LASLlINBI 

16. 

00 

L8XL2INB-1I 

IT. 

00 

LCSL2INBI 

U. 

00 

JH=IEINfil 

19. 

00 

jhP:Jh«i 

20. 

00 

JHM; jh-2 

21. 

00 

00  2S0  NGxl.NGES 

22. 

00 

IGXIGESING.NB  1 

23. 

00 

JGXJGESING.NBI 

24. 

00 

JFrjG-1 

2S. 

DO 

00  210  JiJL.JHP 

26. 

00 

00  210  IXI.MO 

2T. 

00 

210 

Hllfjiro. 

26. 

00 

HII6.JG|:1. 

29. 

DO 

IF  IJF.GT.JHMI  60  TO  226 

30. 

00 

IF  INIYPEUG.JF1.E0. 01  GO  TO  220 

31. 

00 

NDPxNLlNK | NB-1 1 

32. 

00 

00  214  NXI.nOP 

33. 

00 

m=link in.nb-i 1 

34. 

OD 

IF  IICHKIM.NB-ll.EO.IG.ANO.JCHKIM.NB-n.EO.JFt  60  ID  21S 

3S. 

CO 

2  14 

CONTINUE 

36. 

00 

215 

00  21B  NXI.NOP 

37. 

00 

l XL  INK  1 N.NB-I t 

36. 

00 

2  IB 

Ml ICHKIL.NB-1 1 .JF l-RINVl IlB < • N-I 1 *NGE SI *M | 

39. 

00 

2  20 

DO  22S  JXJF.JMM 

40. 

00 

ILtlLFTl J42I-I 

41. 

00 

IR=IRGTIJ«2l-I 

42. 

00 

2  22 

00  22S  I=IL.IR 

43. 

00 

2  25 

Mll«l ,J«2):-IAXtI,J>4HII.J»l>'AYII,JI4HII»l, Jl» 

44. 

00 

A  BBII,J|4HII*I,J*II*CXII.J|4H 11*2, J« 111 

4S  • 

00 

2  28 

00  230  NCxl.NGES 

46. 

00 

I=ICMKINC,NBI-1 

47. 

00 

JXJCHKINC.NBI-l 

48. 

00 

2  30 

rinvila«inc-ii*nges«ngi:axii.ji*hii,j»i i.ayii, ji*hii»i .ji»bbi 

49. 

CO 

A  HII«1,  J«I  I«CXI1  |JI*HII*2.J*1  1 

SO. 

00 

IF  INB.EO.NBLKI  GO  TO  2S0 

SI. 

00 

NOPXNLlNKINBI 

S2. 

DO 

JxIE 1 NB 1 

S3. 

00 

00  240  NXl.NOP 

S4. 

OD 

MZL1NK IN.NB I 

SS. 

DO 

2  40 

RINVIL0*IN-ll4NGES*NGirHIICHK|M,NBI.JI 

S6  . 

CO 

2  50 

HIIG.JGIID. 

57. 

00 

CALL  HATINVlRINVILA.il, N6E S , N6E S 1 

SB. 

CO 

IF  INB.EO.NBLKI  RETURN 

59. 

00 

DO  260  lxl.NGES 

60. 

00 

OD  260  J=1»N0P 

61. 

DO 

R1NV1ILC»IJ-1 I4NGES.II :0. 

62. 

00 

00  26D  KTI.NGES 

63. 

00 

2  60 

R1NV1«LC«<J-1I*NGES.1  IXRINVIILCI  J-ll4NGES.il - 

64. 

00 

A  RINVILA‘IK-1 ) 4NGES.I I4RINVIL0* IJ-1 |4NGES«KI 

6S. 

00 

C  RINV1I1, J.nBI  IS  THE  ‘'ALMOST  HOMOGENEOUS*  *  SOLUTION  *T  THE  J-TH 

66 1 

00 

C  RESIDUAL  POSITION  FORCED  BY  A  RESIDUAL  VALUE  OF  1  AT  THE  1-TM 

67. 

DO 

C  RESIDUAL  POSITION.  HOMOGENEOUS  B.C.’S  ARE  ASSUMED  EVERYMMERE . 

68* 

CO 

C  INCLUDING  THE  TOP  OF  THE  PRESENT  SUBREGION  NB. 

69. 

OD 

JLXJH 

70. 

00 

GD  TO  10D 

71. 

CO 

F  NO 
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VERSION  4 
(continued) 


r 


1. 

01 

SUBROUTINE  NOR (AX, XY.BB, CX ,RINV,R1NV1,0UM0,0UH1,0UM2,F,H,X, 

2. 

00 

11LFT,IR6T,I6CS,JGCS,1Cmk,  JCHK ,NCV .LINK ,HL1NK ,1C ,M0,N0,M2, 

3. 

00 

*  N2,NBLK,NBLK1,MX6S,L1,L2) 

9. 

00 

DIMENSION  iX'H?,N2>,*Y'H2,N2>,BB'H7,N7',CXIH2,N2>, 

S. 

00 

t  rinviii.rinviiii.himo.no), 

(• 

00 

B  Ilf TINOI , 1RGTIN0I ,Ll (NBLK I ,121 NBLK 1 I , 1 GC St HXGS ,NBLK « , 

7. 

00 

C  JGCSIMXGS.NBLK  1  ,1CHK IMxGS.NBLK  I  ,  JCHK 1  MXGS.NBLK  I  ,  NC  V 1  NBLN  1  , 

1. 

00 

0  Ll NK | M2, NBLK1I.NL INK  1 NBLK 1 ) , 1 E 1 NBLK ) 

«. 

00 

DIMENSION  QUMO  1  M2, NBLK 1 I.OUHl 1 MXG S > ,0UH2 I MXGS > ,f IH2.N2) .XI  MO, NO) 

10. 

00 

JSTI 

11. 

00 

00  ISO  NBH.NBLK 

12. 

00 

LC=L2*NB> 

13. 

00 

JF  :lCINB>-2 

1*. 

00 

00  105  J:JS,JF 

IS. 

00 

liniFTl  J«2)-l 

16. 

00 

IRMRGTl  J«7)-l 

17. 

00 

00  105  1-1L, !R 

IB. 

00 

1C5 

Xll«l,J.2):Fll,J)-*xtI.J)*XII,J.ll-*YII,J|9xil«l,JI-BBtl,J)9 

17. 

00 

a  xii*i,j*i)-cxii,j)*xii«2,j*i> 

20. 

00 

IF  IN8.CC.NBLK)  Go  TO  150 

21. 

00 

NGCS=NCVINB ) 

22. 

00 

00  115  NM.NGCS 

23. 

00 

1:ICHKIN,NB>*1 

2*. 

00 

J=JCHKIN,NM)*1 

25. 

00 

I  15 

Ouhknufi  i,  j)-»x(l,J)»x(l  ,j«i>-*yii  ,j>6Xii«i,  J)-BBil,j)» 

26. 

00 

t  Xll«l,J«ll>CXIl,J)9Xt  1*2, J«l)~Xt 1*1, j»2) 

27. 

00 

NOPXNltNKINB) 

26. 

00 

J:1C 1  MB ) 

29. 

00 

00  120  NM.NOP 

30. 

00 

0UM2IN):0. 

31. 

00 

00  11B  M=1,NGC5 

32. 

00 

I  18 

0UM2(N)=0UH2IN)«0UN1 (N)*R1NV1«LC* IN-1 )*NGC5«M) 

33. 

CO 

K:LINKIN,NB) 

36. 

00 

I -ICHK 1 M , NB ) 

35. 

00 

0U«0IN.NB)=X«1, J) 

36. 

00 

120 

x i 1 , j) rx ( l , j )-0uM2 IH) 

37. 

00 

150 

J5- 1C  1 NB ) 

3a. 

00 

00  300  NB 1=1, NBLK 

39. 

00 

NBTNBLK*N8 1*1 

*0. 

00 

JS=1 

91. 

00 

IF  (NB.NE.I)  JS: 1C INB-1 I 

92. 

00 

JF:IC)NB)-2 

93. 

00 

LArillNB) 

99. 

00 

LB:l2iNB-1 ) 

95. 

00 

NGC  5-NE  V  1  NB  ) 

96. 

00 

IF  INB.CQ.NBLK)  GO  10  201 

97. 

00 

j:lC<NB) 

96. 

CO 

NOP  =  NL  INK  1  NB  ) 

99. 

00 

00  200  N-l , NOP 

50. 

00 

NSL INK (N.NB ) 

51. 

00 

1-ICHK tM.NB ) 

52. 

00 

2  to 

XII, J):OUMOlN,NB) 

53. 

00 

2  Cl 

CONTINUE 

59. 

00 

NrlCINB) 

55. 

00 

00  2<T2  j: J5,N 

56. 

00 

00  202  in, MO 

r 

*1 

■1 


rJ 

•i 


i 


H 


»" 


i 

i 

i 


i 

{ 


I 


9 

i 


9 


H 
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VERSION  A 
(continued) 


ST. 

00 

2  02 

Hll.JlrO. 

St. 

00 

00  210  NII.NGES 

59. 

00 

I=ICHMN,NB»-1 

60. 

00 

j:jchhin,nbi-I 

61. 

00 

210 

OUHl  mi:Fl  I  ,J)-AXI  I  ,  J)9XI  I  ,  J»I  )-AY  (  I  ,J)9X  1  I.  1  ,  J)-B8I  I,  J|9 

62. 

00 

A  XI1*1,J*I »-CX<I,JI*XII»2, J.lt-XI I«I,J*2> 

61. 

00 

00  220  Ntl  ,nGES 

6*. 

00 

0UH2IN>=0. 

66. 

00 

00  218  h:1,NGES 

66. 

00 

211 

ouh2(nhouh2(ni.ouhhhi»rinv  ila.in-i  i*nges«h> 

67. 

00 

I=lGESIN,NBt 

66. 

00 

JXJGESIN.NBI 

69. 

00 

Htl, J|:0UH2(N» 

70. 

00 

2  20 

XtI,J)=XII,Jt*0uM2«Nl 

71. 

00 

IF  INB.EO.II  GO  TO  250 

72. 

00 

00  222  n:1,HX6S 

IS. 

00 

2  22 

0UH1 I N)-0. 

79. 

oo 

NOPRNLlNM  N8-I  I 

75. 

oo 

J:JS 

76. 

00 

00  230  N-I .NOP 

77. 

00 

h:linkin,nb>H 

78. 

00 

I:ICHKIH,NB-I 1 

79. 

oo 

2  70 

OUHIIHUHd,  J»lt 

80. 

oo 

C  WH£N  80UN0ARV  IS  IRREGULAR,  THERE  IS  SOHE  WASTEO  CALCULATION  1 

81. 

oo 

C  OROER  TO  AVL10  LOGICAL  0CC1S10NS  OR  EXTRA  STORAGE  INGES  COuLO 

82. 

00 

C  REPLACEO  BY  NLINKINB)) 

81. 

00 

NGESXNEVINB-I 1 

89. 

00 

Oo  290  N=I,NOP 

85. 

00 

OuH2(Ni:0. 

86. 

oo 

00  238  YCI.NGES 

87. 

00 

2  38 

0UH2I N>rou*2IN».0u»I IN  HRINVI  1  LB*  IN-I  XNGES.HI 

88. 

00 

«=LINK 1 N.NP-I ! 

89. 

oo 

2  90 

MIICHK(HtNB-Il ,JCHX (H.NB-I I >:0UH2INI 

90. 

00 

2*0 

00  300  JRJS, JF 

91. 

oo 

IL=ILFT(J.2|-I 

92. 

oo 

IRI1RGT  t J«?1«I 

93. 

00 

00  300  im.IR 

99. 

00 

H«I«I,J.2l--(AXIItJ|9HII,J.l».AY(I,J19HII.l,J»» 

95. 

00 

A  BB(I,JI9HII*I,J.1}.CX(Z,JI9H(1.2,J.I)) 

96. 

00 

ICO 

X( 1*1, j.?irx(IM,J*2l*HiI«I, J»2I 

97. 

oo 

ENO 
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VERSION  A 
(continued) 


I. 

01 

SUBROUTINE  TOPOL • 1LE T, TR6T, I 6ES.J6ES, I CHK, JCHK, NEV, LINK, 

2. 

00 

1NLINK,NTVRE,IE,H,N,H1 ,NI,H2,NBIK,NBlk1 .HXGS.LI ,L2,L01,L02> 

3. 

00 

01  hen  SI  ON  ItrilNt ,IRGTIN),I6ESIHX6S,NBLK>, 

9. 

00 

A  JGESIHAGS.NBLKI.lCHKtHXGS.NBLK  >, JCHK 1 HAGS , NBLK I , NE V 1 NBLK t , 

5. 

00 

BL INK  IH2.NBLK1 ) .NLlNKlNBLKIt .NITRE  1 H ,N > , I E I NBLK • , L I I NBLK ) , L2I NBLKI > 

t. 

00 

00  10  J=2,NI 

7. 

00 

1=0 

8. 

00 

2 

1=1*1 

9. 

00 

ir  INTYPEtl«JI*NTyR£ll,J-I).EO.O.ANO.I.LT.HI)  GO  TO  2 

10. 

00 

iLrTtj»=i 

II. 

00 

1  =H 

12. 

00 

9 

1=1-1 

13. 

00 

ir  INTYPCIl,J)*NTYPEtl,J-I>.Ee.0.AN0.1.6T.l>  GO  TO  9 

19. 

00 

1C 

1RGT I J • =1 

IS. 

00 

IRGTt 11=2 

1  b  • 

00 

iLrrti»=2 

17. 

00 

iLri<N>=2 

IS. 

oo 

IRGTIN>=2 

19. 

00 

00  IS  J=2.N1 

20. 

00 

ir  IlLETI jl.LE.IRGTtJI)  GO  TO  IS 

21. 

oo 

IRGT(J)=2 

22. 

00 

iLrnj»=2 

23. 

00 

IS 

CONTINUE 

29. 

00 

JHIN=2 

25. 

00 

HALNK=0 

2b. 

00 

LI  1 1 1=0 

27. 

oo 

L2I 11=0 

28. 

00 

00  100  NB=I,NBLK 

29. 

00 

NGES=0 

30. 

00 

NCHK=0 

31. 

00 

JHAX=1£<NB> 

32. 

00 

00  90  1=2, HI 

33. 

00 

NH=NTYRE II , JHIN-I I 

39. 

00 

NU=NTYREII.JH1NI 

35. 

00 

00  90  J=JHIN.JHAA 

3b. 

00 

NL  =  NH 

37. 

00 

NH=NU 

38. 

00 

NU=NTVPCI1,J«1» 

39. 

00 

ir  1 NH .EO .91  GO  TO  90 

90. 

00 

ir  inl.eo.i.ino.j.ne.jhini  go  to  20 

91. 

00 

N6ES=NGES* 1 

92. 

00 

ir  inges.gt.hagsi  go  to  200 

93. 

00 

16E SINGES. NB 1=1 

99. 

oo 

JGE  SINGES, NBI  =  J 

95. 

00 

2  C 

ir  INU.EO.I.ANO. J.NE.JHAX)  GO  TO  90 

9b. 

00 

NCHK=NCHK«1 

97. 

oo 

ICHk INCHK ,NBI=1 

98. 

00 

JCHKI NCHK.NB l=J 

99. 

00 

9  C 

CONTINUE 

SO. 

00 

NEWINBI=NGC$ 

51. 

00 

ir  INB.EO.NBLKI  GO  TO  100 

52. 

00 

NLNK  =0 

53. 

00 

00  BO  K=1,NGES 

59. 

00 

1=ICHKIK,NB> 

55. 

00 

J= JCHK IK ,NB • 

5b. 

00 

ir  INTYPEIl,J«ll.CO.O>  GO  TO  80 

57. 

00 

NLNK=NLNK«I 

5a. 

00 

LINK  1 NLNK ,NB l=K 

59. 

00 

ac 

CONTINUE 

bO. 

00 

NLINKINB  »=NLNK 

bl. 

oo 

nalnk=haaoihalnk ,NLNK*NGES) 

82. 

00 

LIINB*It=LllNBI*NGES**2 

83. 

00 

ir  INB.NE.NBLKlt  L2|NB‘1»=L2INB»*NGES*NLNK 

89. 

00 

LRINWI=L2INB»«NGES*NLNK 

85. 

00 

ICO 

JH1N  = JHA A* 1 

bb. 

00 

c  If 

halnk,gt.nges*nges.  there  is  sohe  wasteo  storage  that  is  vert 

87. 

00 

c  oirricuii  to  prograh  arouno,  because  or  the  dual  role  or  the  list 

88. 

00 

C  BLOCK  or  R I NV  ,  WHICH  15  USED  IS  8  SCRATCH  ARE*  EOR  RINVI  OATA. 

89. 

00 

lrinv=liinblkmhaa9ihalnk,nges9*2> 

70. 

00 

WRITEIb.159)  LRINV.LRINVl 

71. 

00 

ISO 

rORHAT I55H  01HENSI0NS  REOUIREO  TOR  RINV  ANO  RINVI  ARE  L01,  L02  =  . 

72. 

00 

AI5.2A.15.SH./129H  ir  EITHER  INPUT  VALUE  LOI  OR  LOT  IS  TOO  SHALL.  T 

73. 

00 

BHE  COHPuTATION  IS  HALTEO  HERE,  ANO  THE  USER  HUST  INCREASE  LOI  ANO/ 

79. 

00 

COR  LOT  T0/72H  THE  INDICATEO  VALUES.! 

75. 

00 

ir  IIlOI.LT, LR1NVI. OR. ILOT.LT.LRINVIM  STOP 

7b. 

00 

RETURN 

77. 

09 

2  CO 

WRITE IG.202T 

78. 

00 

2E2 

roRHAT 1 92h  0IHENS10N  PARAHETER  hags  IS  TOO  shall  rOR  SPECiriEO  TOP 

79. 

00 

AOGRAPHY  ANO  SUBREGION  SPEC 1F1CA T10N. * 

80. 

00 

STOP 

81. 

00 

eno 
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VERSION  4 
(continued) 


00 

SUBROUTINE  M»TINV(B,N,HI 

00 

OIMENSION  BIN,11iB1I100*iB2I100> 

00 

H1=H-1 

00 

00  110  1=1. Ml 

00 

Bl(ll=l«/B(l«ll 

00 

BIT. 11=1.0 

00 

00  112  J=1.N 

00 

112 

B(I,JI=B(I,J1*B1«11 

00 

1P1=1*1 

00 

00  120  I1=1P1.N 

00 

120 

B1(I11=B(I1.11 

00 

DO  125  I1MP1.M 

00 

1  25 

BUI. 11=0. 

00 

DO  127  J=1.H 

00 

127 

B2I J1=BI I , Jt 

00 

DO  135  11=IP1.N 

00 

DO  135  J=1 . M 

00 

115 

B(I1,JI=B(I1,JI-B1(I11*B2(JI 

00 

1 10 

CONTINUE 

00 

81(11=1. /B(H, HI 

00 

6 IH.H 1=1 . 

00 

00  1*0  J=1,H 

00 

1«0 

BIH,JI=6(H,JI*B1I1I 

DO 

DO  150  1=2, H 

00 

DO  155  12=1,1 

00 

1  55 

B1II21=B<12,1 1 

00 

IHI  =  I-1 

00 

00  156  12=1, INI 

00 

1*6 

6(12.11=0. 

00 

00  157  J=1,H 

00 

1  57 

821 J1=BII,JI 

00 

181=1-1 

00 

00  160  12=1, TUI 

00 

DO  160  J=1,H 

00 

1  60 

B 1 12, J1=BI 12. Jl-Bl (121*621 J> 

00 

1*0 

CONTINUE 

00 

ENO 

40 


VERSION  5 


1* 

02 

PR06RAN  JAYNOR  1 1NPU1  ,  OUTPUT,  TAPESUNPuT  .TAPE6I0UTPUT  1 

2. 

00 

01  pension  AXI7,9I,ATI7,9I,OOI7,9I,CXI7,9I,OI7,9I,FI7,9I, 

3* 

00 

A  RTNVI7,7,3| ,R1NV1I7,7,2I,0UN0I7,2I,0UH1I7I,0UH2I7I,XI9,11I, 

l« 

00 

0  HI9,11I,ERRI9,11I,XXI9,11I,1EI3I 

S. 

00 

DATA  H,N,NBLK/9, 11 ,3/ 

6  » 

00 

DATA  IE/9,7, 10/ 

t. 

00 

C  OI£TRTCH*S  N001 F 1 E  0  HAO«LA  SOLVER .  VERSION  5. 

o. 

00 

(  THIS  VERSION  15  IDENTICAL  TO  VERSION  «,  EXCEPT 

*. 

00 

C  IRREGULAR  BOUNDARIES  ARE  NOT  ALLOWED.  IT  IS  ALSO  EOUTVALENT  TO 

10. 

00 

C  HAOALA'S  BASIC  SOLVER.  ALLOWING  NON “SEPARABLE  OPERATORS  ON  A 

11. 

00 

c  non-uniforh  he sh  with  regular  boundaries. 

12. 

00 

C  THE  E0UAT10N  TO  BE  SOLVED  IS 

13. 

00 

C  AXll.JI  9  XII,  J.  11  «  CXU.JI  *  X  tl«2,J*l  I  *  BBll.JI  9  XI1«1,J«1I  ♦ 

1*. 

00 

C  AVIl.JI  9  XI1M.JI  ♦  XtI«l,J»2l  =  F 11, Jl, 

IS. 

00 

c 

111=1, n-2I,Ij=1,n-2Ii 

16. 

00 

C  WHERE  XII, Jl,  XIH.JI,  XII, 11,  AND  XII, Nl  ARC  LEFT  UNCHANGED. 

IT. 

00 

H1=H-1 

16. 

30 

ni=n-i 

19. 

00 

H2IH-2 

20. 

00 

N2=N-2 

21. 

00 

NOLXIXNOLX-I 

22. 

00 

00  100  J=1,N2 

23. 

00 

00  100  1=1, H2 

2*. 

00 

AX ll,J 1=1. 9.019 11 -91992 

25. 

00 

CXI1,J|=1.-.01*|1-9I**2 

26. 

00 

AT  11  ,JI  =  1.9.01 9| J-51 992 

27. 

00 

OBI I,JI=I-AXI1 ,JI-AV!l,JI-ex 11, Jl 9.01*1 J-S 1992-1. 11/ 

26. 

00 

A  ll«-.01*IJ-SI**2l 

29. 

00 

100 

continue 

30. 

00 

0£ NON: I.2S9H9N l**2 

31. 

00 

00  170  J  =  1 ,N 

32. 

00 

00  170  1  =  1, H 

33. 

00 

X II, Jl=ll-1. I911-R19IJ-1. 1*1 J-NI/DENOH 

3*. 

00 

XXII, Jl=0. 

35. 

00 

IF  iI.EO.l.OR.l.EQ.M.OR.J.EO.l.OR.J.EQ.NI  XX 1 1 , J l=x 11 , J I 

36. 

00 

170 

CONTINUE 

37. 

00 

FSUH=0. 

36. 

00 

00  171  J=2,N1 

39. 

00 

ISJ-l 

*0. 

00 

00  171  1=2,H1 

*1. 

00 

X  =  l-I 

*2. 

00 

fix,li:ixik,li*xii,j>»avik,li*xii,ii»cxik,li*xii»i.ji* 

93. 

00 

A  XII,J«ll«BBIK,LI*XlI,JI 

•9. 

00 

1  71 

F  SUH  =  F SUH.A  BS 1 F 1 R ,L 1 1 

95. 

00 

FSUH=H2*N2/FSUH 

96. 

00 

TO-SECONOlTIHI 

97. 

02 

CALL  JATIAX,AT,BB,CX,R1NV,R1NV1,H,TE,H,N,H2,N2,NBLR,NBLK1I 

9*. 

00 

TIHSUHIO. 

99. 

00 

00  101  J=2,NI 

53. 

00 

00  101  1=2, Nl 

51. 

00 

011-1, J-1I=FI1-1,J-1I 

52. 

00 

1  01 

XIT.JISO, 

53. 

00 

AllNEsO. 

59. 

00 

T0=SEC0N0ITJHI 

55. 

02 

CALL  NOR IAI,AT,BB,CX,R1NV,RZNV1 , OUNO • OUN 1 , 0UN2 ,0,N,XX,1E,N,N,N2 

56. 

00 

A  N?,NBLK,NBLK1I 

57. 

00 

00  990  J=2.NZ 

56. 

00 

00  990  1=1, Nl 

59. 

CO 

9  50 

XII,JI=XII,JI'XXI1,J1 

60. 

00 

00  991  J= 1 , N 

61. 

00 

00  991  1=1, N 

62. 

00 

951 

xxii,  ji=n. 

63. 

00 

00  992  J=2,N1 

69. 

00 

l=J-l 

65. 

00 

00  992  1=2, HI 

66. 

00 

R  =  I-1 

67. 

00 

952 

0 IK,LI=riK,LI -AT IR,L 1*XI1, J-ll-AXIR.LI 9XII-1, Jl -00  IK ,L 19X11, Jl- 

6t. 

00 

A  CXIK,ll*Xll«],JI-Xll,J«ll 

69. 

00 

iroo 

CONUNUE 

70. 

00 

9  CS 

FORHATI1X,1P15E0.1I 

71. 

00 

ERSUNSO. 

72. 

00 

00  790  1=2, HI 

73. 

00 

R  =  l-1 

79. 

00 

00  790  J=2,N1 

75. 

00 

L  =  J-I 

76. 

00 

ERR  11 ,JI=FIR,L l-ATIR.L 19X11, J-ll-AXlK.l 19X11-1, JI-BBIR.ll* 

77. 

00 

A  XI1,JI-CXIX,LI*XI1«1,JI-XI1,J«1I 

7t. 

00 

ERSun=ERSuN«ABS|ERRI1,JII 

79. 

00 

7  90 

CRRIJ,JI=EPRll,JIRFSUH 

•  0. 

00 

rRSUN;ERSUH,FSUN/|Nr*N2l 

oi. 

00 

WRITE  1 6,909 1 

82. 

00 

9C9 

FORHATI/9CV,?OHNORNAL17EO  RESIDUALS  1 

03. 

00 

00  792  1=2, HI 

(9. 

00 

WRITE  16,9051 IERRII,JI,J=2,N1I 

•  5. 

00 

792 

CONTINUE 

06  . 

00 

WRITE! 6,906 1  ERSUH 

07. 

00 

9E6 

F0PHA1I/3TH  NEAN  ABSOLUTE  NO»NA112EO  RESIDUAL  =  , 1PE9.2 1 

00. 

00 

ENO 
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VERSION  5 
(continued) 


SUBROUTINE  JAVIAX,AT,BB,CX,R1NV,R1NV1,M,1E  ,*0  ,N0  .HZ,  N2  ,NBL*  , 
1NBLKII 

DIMENSION  AX IM2,N2> ,AYIM2,N2> ,BBf H2.N2) ,CXIM2,N2I ,R1NVIM2,H2,NBLK> 
A,RINV1IM2,N2,NBLKI  I  .H(M3,nDI  .ieinblm 
JL  =  1 
NB-0 

ICO  NB=NB» 1 

JH=1EINB> 

JHP=JH«1 

JHM;jM-2 

JG=JL«I 

00  2S0  NG=1,H2 
1G=NS«I 

00  210  J=JL, JHP 
00  210  1=1. HO 
210  Htl,J):0. 

HI  IG  .  JG 1=1. 

IE  INB.EO.II  60  TO  220 
00  21B  N=1,H2 

218  HIN«I,JL)=R1NV1ING,N,NB->1> 

220  00  225  J= JL . JHM 

00  225  1=1, M2 

2  25  HII«I,J«2I=-IAX|I,JI»HI1,J»1>«ATI1,J>*H|1«1,  JI*BBII.JI*Htl«llJ«ll* 

A  CXtl,JI*Hll«2,J«l>> 

J=JH-1 

00  230  1=1, M2 

2  30  R1NV(NG,1,NB I =A  X  1 1 , J I *H II ,  J*  I I  ♦  AY (I ,  J 1 4  H I 1*1, Jl *BB  I 1 ,  J I  * 
k  HII«1, J*1I«CXI1,JI*H|I«2, J»I I 
IE  INB.EO.NBLKt  GO  10  2SD 
J=IE INB I 
00  240  N=1 , M2 

240  RINVING.N.NBLK l=HIN»I,Jt 
250  H I 1G , JG ) =0 • 

CALL  MAIINVIRINVI1,I,NBI,M2,M2I 
IE  INB.EO.NBLKI  RETURN 
00  260  1=1, M2 
00  260  J=1 ,M2 
R1NVII1, J,NB>=0. 

00  260  K=1,M2 

260  R1NV1I1,J,NBI=R1NV1II,J,NB>-RINVII,K,NBI*RINVIK,J,NBLK> 

C  RINVIII.J.NB)  15  THE  ''ALMOST  HOMOGENEOUS"  SOLUTION  AT  THE  J-TH  OPEN 
C  RE510UAL  POSITION  FORCED  BT  A  RESIDUAL  VALUE  OF  1  AT  TmE  I-TH 
C  RESIDUAL  POSITION.  HOMOGENEOUS  B.C.'S  are  ASSUMED  EVERTHHERE, 

C  INCLUDING  TME  IOP  OF  THE  PRESENT  5UBREG10N  NB . 
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